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We consider a general five-dimensional sigma-model coupled to gravity, with any number of scalars 
and general sigma-model metric and potential. We discuss in detail the problem of the boundary 
conditions for the scalar fluctuations, in the case where the fifth dimension is compact, and provide 
a simple (and very general) algorithmic procedure for computing the spectrum of physical scalar 
fluctuations of the fully back-reacted system. Focusing in particular on the conditions under which 
the spectrum of scalar excitations (glueballs) contains parametrically light states, we apply the 
formalism to some especially simple toy models, which can be thought of as the gauge/gravity duals 
of strongly-coupled, non-conformal four-dimensional gauge theories. Our examples are chosen both 
within the context of phenomenological effective field theory constructions (bottom- up approach), 
and within the context of consistent truncations of ten-dimensional string theories in the supergravity 
limit (top-down approach). In one of the examples, a light dilaton is present in the spectrum in 
spite of the presence of a bad naked singularity in the deep IR, near which the RG flow of the dual 
theory is certainly very far away from any fixed point. If this feature were to persist in a complete 
model in which the singularity is resolved, this would prove that a light dilaton is to be expected 
in at least certain walking technicolor theories. We provide here all the technical details for testing 
this statement, once such a complete model is identified. 
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I. INTRODUCTION 



Many physical systems are described by strongly-coupled field theories, the dynamics of which is encoded in the 
fixed points of their renormalization group flow. Long-distance properties of such systems can be classified in terms 
of universal coefficients, which depend only on the properties of the system very close to the fixed points, but do not 
depend on the model-dependent features of the flows away from such fixed points. On the other hand, there are many 
cases in which the knowledge of the physics at the fixed points does not provide enough information as to allow one to 
compute phenomenologically important physical quantities that are experimentally measurable. In these latter cases, 
traditional field theory techniques are not powerful enough to yield robust predictions, which can be compared to the 
experimental data, mainly because of the strongly-coupled nature of the underlying dynamics. 

One such example emerges in the context of dynamical electro-weak symmetry breaking, in particular in walking 
technicolor [l[. In this case, the underlying strong dynamics is quasi-conformal (approaching an IR fixed-point) over a 
range of energies above the electro-weak scale, but ultimately yields to confinement and to the formation of symmetry- 
breaking condensates at the electro-weak scale. One might expect that there exists a sense in which the condensates 
break spontaneously the (approximate) dilatation symmetry of the system near the (approximate) fixed point, thus 
leading to the appearance of a light dilaton (the pseudo-Goldstone boson of scale-invariance) in the spectrum of 
composite states. From a phenomenological point of view, this example is of most urgent importance, because such a 
light dilaton might mimic the properties of the Higgs particle of weakly-coupled models such as the minimal version 
of the Standard Model 0. Unfortunately, because of the strong dynamics, and because the physics of a massive state 
such as the dilaton depends not only on the (universal) properties of the fixed-point, but also on the (non- universal) 
RG flow that yields confinement and chiral symmetry breaking itself, neither a firm confirmation nor a disproof of the 
existence of a light dilaton in walking technicolor has been provided by conventional field-theory methods, in spite of 
many attempts [3|. For recent work supporting the idea that such a light dilaton exists, see for instance [j, @. 

In recent years, the discovery of gauge/gravity dualities provided a new tool, that allows to reformulate field-theory 
problems emerging within four-dimensional strongly-coupled systems in terms of weakly-coupled extra-dimensional 
systems [||. In its original formulation Q, the idea is to relate a particularly simple and symmetric 10-dimensional 
background (AdS$ x S 5 ) to a very special conformal four-dimensional theory (Af = 4 supersymmetric SU(N) gauge 
theory) , in the sense that a prescription is given for computing the generating functionals of correlation functions on 
the two sides of the correspondence, and the physical results agree. More recent developments provided large classes 
of dual models that correspond to non-conformal field theories with much less supersymmetry, such as those in which 
the 10-dimensional background is constructed starting from the conifold and its variations [8j, and those that are 
related to controllable deformations of the Af = 4 field theory [9j . 

It turns out that, for several reasons, all the models that are of phenomenological interest share some very general 
properties. In particular, the 10-dimensional metric is always written in terms of a non-compact five-dimensional part 
(four directions of which are directly related to the dual four-dimensional space, with the fifth dimension related to 
the energy scale at which the dual theory is tested), and a compact (internal) five-dimensional space, the isometries of 
which are related to the internal global symmetries of the field theory. Formally, this means that it is often possible, 
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and very convenient, to study the 10-dimensional dynamics by first studying its five-dimensional reduction. One also 
implements a consistent truncation that reduces the number (and simplifies the action) of the resulting active fields. 

The problem is thus reduced to finding a background solution of the classical equations of the truncated five- 
dimensional theory, and then use it to construct the lift to a background solution of the full 10-dimensional system. 
In practice, this means that what one has to solve, in order to determine the background, are the classical equations of a 
five-dimensional sigma-model of n scalars coupled to gravity. Since one usually wants to preserve the Lorentz structure 
of the four-dimensional space-time, one also assumes that the background depends only on the radial direction. With 
all of this, the complexity of the original problem of finding fully back-reacted 10-dimensional backgrounds is turned 
into the more treatable problem of solving a one-dimensional classical system. 

One can even use this formalism dispensing with the original problem of finding a consistent truncation, by sim- 
ply writing a five-dimensional phenomenological model that captures the most important features of the dynamics, 
postponing the problem of its completion (all we are going to say applies also in the fake supergravity context [lolp. 
In this spirit, a vast literature of phenomenological models exists which aim at capturing the most important aspects 
of the gauge/gravity dualities without dealing with the technical difficulties of a complete string-theory construction. 
Relevant examples include the Randall-Sundrum model [ll|, in which the sigma-model consists just of a cosmolog- 
ical constant, the Goldberger-Wise mechanism [l2[ (see also fl3T0 . in which the model contains only one scalar (see 
also fl4|), a nd many applications, such as the Higgsless models |15| . the AdS/QCD models [lfj], some composite- Higgs 
models [TtJ and the holographic technicolor models [l8| . Recently, a similar strategy has been proposed also in order 
to study lower-dimensional condensed matter systems (see [l9| for an introduction to the subject). 

After reinstating (in the five-dimensional action) the dependence on the Minkowski directions, one can study the 
spectrum of classical fluctuations, which can be interpreted as the composite states of the dual theory (the glueballs, 
for instance). From this, one can finally access those very non-trivial properties of the strong dynamics that we started 
by describing in the beginning of this introduction, and ask hard questions such as whether a specific model yields a 
light dilaton in the spectrum. Yet, there are still two difficulties to overcome, before a final answer to these questions 
can be provided. 

First of all, because the spectrum of massive states is not a universal property, one has to construct and study 
a variety of explicit models, possibly such that a lift to a complete 10-dimensional theory exists, and such that the 
dual field theory has all the properties required by phenomenology. In the specific context of walking technicolor, 
this program has recently been initiated, with some very encouraging results [H, 120I l2~lj (see also [22j for a more 
extensive discussion of what one might want to achieve along this line). Yet, at present we are nowhere close to 
having constructed the actual 10-dimcnsional dual of a phenomenologically relevant four-dimensional model, as a 
substantial amount of model-building is required in order to do so. 

A second difficulty is of a more technical nature, and is the main subject of the present paper. When studyin g th e 
spectrum and the properties of the fluctuations, a long preliminary work appears to be necessary (see for example [231 ]) 
— mainly because of the non-trivial mixing between fluctuations of the bulk scalars and the five-dimensional metric. 
One hence would need a formalism that is general and simple enough to correctly incorporate the relevant dynamics 
without having to analyze all the details on a model-by-model basis. In the case of a single scalar with trivial sigma 
model, this program has been addressed time ago by several collaborations (see for example [HI)- Yet, the formalism 
developed by these authors needs to be extended far beyond the level needed for a single scalar, in order for it to 
apply to a realistic gravity dual of a strong dynamics. 

In pBj . a step towards a systematic resolution of this technical problem was taken. By introducing appropriate 
gauge- invariant combinations of the original fluctuations (see also [26j and [HI), which from the four-dimensional point 
of view correspond to scalars, vectors and tensors, it was shown that, given a completely general (two-derivatives) 
sigma-model with n scalars, and with a superpotential W (from which the 5d potential V can be derived), it is possible 
to algebraically manipulate the system of linear differential equations so as to rewrite it in a sigma-model covariant 
form and reduce it to a set of n second-order equations for the same number of physical fluctuations, from which the 
spectrum of the spin-0 sector of the theory can be derived. It was also shown in [27| that the formalism can easily 
be generalized to the case where the superpotential W is not known (or does not exist), in which case one needs to 
know the sigma-model metric and the potential. 

This important formal tool works thanks to the fact that one can use the five-dimensional diffeomorphism invariance 
in order to remove some of the unphysical fluctuations. However, in practical applications one has to generalize this 
instrument further, so that it applies to the case where the radial direction is not infinite and, hence, boundary 
actions may need to be added. For example, in many cases an IR boundary is present because of an end-of-space 
in the geometry (which must be the case when discussing the dual of a confining gauge theory). Also, in the UV it 
is often necessary to work with a finite cutoff, for three possible generic reasons. It might be known that the dual 
field theory requires a UV completion above a given scale, so that the UV cutoff is actually physical. Retaining 
a UV cutoff is also necessary for technical reasons related to holographic renormalization [2g|. And finally, it may 
be that the strongly-interacting dual theory is (weakly) coupled with an external (weakly-coupled) four-dimensional 
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sector, modeled by UV-boundary interactions. In all these three cases, the boundary terms are going to break the 
five-dimensional diffeomorphism invariance, and hence some caution has to be used when applying the gauge-invariant 
formalism. 

In this paper, we provide the general form of the boundary conditions for the scalar fluctuations, both in the 
case in which the five-dimensional dynamics is known in terms of a superpotential, and when only the potential 
exists, without any restriction on the number of sigma-model scalars or on the sigma-model metric. We discuss the 
residual freedom in the form of the boundary conditions, in particular in relation to the light scalars in the spectrum, 
one (linear combination) of which may be interpreted as a dilaton, while any others correspond either to ordinary 
pseudo-Nambu-Goldstone bosons of approximate global symmetries of the dual theory, or are the result of accidental 
cancellations. We illustrate the formalism hence derived by applying it to a set of simple phcnomenological examples, 
for which we provide both (approximate) analytical results for the spectrum, and (exact) numerical studies. 

We do not treat here the problem of holographic renormalization, which we postpone to future work. In particular, 
we will always consider the models as defined on a compact fifth-dimension, away from any possible singularities. 
In this way, all the states are going to be physical. In practical examples, one has also to decide how to couple the 
model to possible external (weakly-coupled) systems, and how to take the limits in which the IR and UV boundaries 
are removed. We will briefly comment on these issues in due time, but we are not going to provide a systematic 
prescription for doing so. 

A. The algorithm 

The main purpose of the paper is to provide the reader with a simple algorithmic procedure for computing the 
spectrum of scalar excitations. Assuming that a five-dimensional model is of interest (irrespective of the fact that it 
is built either as a phenomenological model, or as the consistent truncation of a given supergravity, or superstring, or 
M-theory), one should go through the following steps. 

• One must first ensure that the model can be written in the general form we provide, i. e. as a two-derivative, 
5-dimensional action involving n real scalars coupled to gravity. We do not consider the case where higher- 
derivative terms are present, and ignore the possibility that higher-spin fields (such as gauge bosons or fermions, 
for instance) are relevant in determining the background. 

• One must assume that two boundaries are present in the fifth dimension, representing the UV and IR cutoffs 
of the dual theory. These cutoffs may have a physical meaning, in which case all the results will explicitly 
depend on the dual scales. Or they may be thought of as regulators, in which case the physical results should 
be obtained by extrapolating the final results of the calculation to the actual physical case (typically, one wants 
to remove the UV cutoff completely, while the IR cutoff will approach the end-of-space). The presence of the 
boundaries means that boundary actions must be added. We provide the most general form of such boundary 
actions, subject to the limitation that we do not include terms that depend on the Minkowski four-momentum 
q 2 . As we will see, this form is very constrained, although we restrict ourselves to the quadratic order. 

• One has to solve the system of background equations, and find a suitable background, with the ansatz that 
all the background functions depend only on the radial direction. We provide the general form of the bulk 
equations and their boundary conditions derived from the complete action. We do so both in the case in which 
the sigma-model is described in terms of a potential, but also in the case in which a superpotential description 
is known. The latter has the advantage that the background is completely determined by a set of first-order 
(coupled and non-linear) differential equations. 

• One has to solve the linearized, second-order equations for the fluctuations around the background. The spec- 
trum is determined by solutions that satisfy the boundary conditions both in the UV and in the IR. We provide 
the complete set of bulk equations [25| . directly in the physical basis, written in terms of the background solu- 
tion. And we provide the boundary conditions, the form of which depends again on the background functions, 
but also on a set of parameters that incorporate the residual freedom in the choice of boundary actions. These 
parameters should be chosen on the basis of physical principles, and are hence model-dependent. However, in 
the absence of symmetry reasons, there is only one such choice that ensures the absence of any fine-tuning in 
the physical results, and in this limit the boundary conditions depend, again, only on the functions determining 
the background, evaluated at the boundaries. 

This procedure is very general, and since we write the fluctuation equations and boundary conditions already in 
terms of physical fields, no algebraic manipulations are needed. The reader who wants to apply this procedure can 
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directly write the final system in terms of n gauge-invariant fluctuations and solve it (in most of the cases, this must 
be done numerically). The formalism of [25[ allows to write all the relevant equations in an elegant form that is 
fully covariant with respect not only to the space-time, but also the internal sigma-model geometry. One should be 
careful in correctly using all the covariant derivatives, which are determined both by the space-time and sigma-model 
connections: we will present all the relevant (heavy) notation in Section [TT1 

Again, we must stress that the generality of the boundary terms for the fluctuations is limited by the fact that 
we do not include q 2 -dependent terms. These are model-dependent, and important in the context of holographic 
renormalization, when trying to take the UV cutoff to infinity, a problem that we postpone to future work. 

B. Reader's guide 

The paper is organized as follows. In Section [TT| we summarize the basic formalism we use. This is a rather 
technical section, which may be skipped at first reading. However, all the material contained is necessary in order to 
correctly interpret and use the basic equations appearing in the paper, and we find it convenient to group together 
all the necessary definitions in one place. Also, we discuss here some subtleties emerging in the introduction of 
gauge-invariant variables in the presence of boundaries, which clarify and complete the literature on the subject. 

In Section Hm we write the action of the sigma-model coupled to gravity, the equations determining the background, 
and the final differential equations and boundary conditions satisfied by the physical fluctuations. The derivation of 
these results is summarized in the appendices. All the relevant equations are written in terms of the sigma-model 
metric, the background fields, and the potential (and, when available, superpotential). 

In Section llVI we present three particularly interesting examples, and derive some analytical results. In Section IVl 
we apply the mid-point determinant method [I^] to study these examples numerically. The main purpose of the 
examples is to illustrate the procedure, and hence we choose them to be particularly simple. However, it turns out 
that their physical interpretations are quite interesting, and that this set of exercises also provides some important 
insight into how the regulation procedure may work in non-trivial physical cases. 

The first example is based on the same action used in the GW mechanism, and allows us to compare our results 
to the literature, but also to generalize the results and discuss many interesting subtleties that have been ignored in 
the past. In particular, we explicitly show that some freedom in the definition of the boundary terms results in the 
possible appearance of additional light scalars besides the dilaton, and that hence one has to exercise some caution in 
interpreting the results. We also show that a light dilaton is present (in great generality) not only when the scaling 
dimension A of the field-theory deformation encoded in the background is small (A <§C I), but also for any A > 2. 

The second example is taken from a peculiarly simple five-dimensional model constructed by consistent truncation 
of type-IIB supergravity. The dual gauge theory has many properties that resemble those of a QCD-like theory, in 
the sense that the formation of a condensate in the IR takes the theory away from its fixed-point, presumably leading 
to confinement. Unfortunately, the model (studied here at zero temperature) suffers from the appearance of a naked 
singularity in the background, which limits its physical meaning. Yet, it is interesting to study what happens to the 
spectrum in the limit where the IR cutoff approaches the singularity. As we will see, the procedure adopted here yields 
a spectrum that, while distorted by the presence of the singularity, does not show any signs of pathologies, suggesting 
that the procedure that we follow removes some of the unpleasant features of the background at the singularity. 

The third example is a phenomenological model yielding a background that can be interpreted in terms of the RG 
flow between a UV fixed-point and an IR fixed point. We study in some detail what happens to the spectrum by 
comparing several backgrounds that differ only by the value of the scale at which the transition from the proximity 
to one fixed point to the other takes place. 

In Section IVT1 we present a set of field-theory arguments aimed at explaining the results of the previous two sections. 
We elaborate on possible interpretations of the five-dimensional models in terms of dual, strongly-coupled theories, 
and derive some lessons about the physics of the three models we considered. These lessons extend to any model the 
examples somehow approximate, and are hence of general interest. 

In Section I VIII we conclude, by summarizing the most important equations needed in the proposed algorithmic 
procedure, by commenting on the limitations and subtleties involved in using the algorithm itself, by summarizing 
briefly the physics lessons we learned, and finally by outlining some possible future applications. 

II. FORMALISM 

We introduce here the main definitions and conventions we use in the paper. We do not make explicit use of 
supergravity transformation and other supersymmetric properties. We start from the definition of the geometric 
properties of a sigma-model of n scalars coupled to gravity in five dimensions. Most of the notation and conventions 
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we use are taken from (25j . The notation is somewhat heavy, and hence we devote some time to explain it, explicitly 
showing all the definitions used in the whole paper. We then explicitly discuss the effect of gauge transformations 
and introduce the gauge invariant variables that will be used throughout the paper. 



A. Geometry 

All the equations we write make use of the geometric properties of the sigma-model, and are hence completely 
covariant. The space-time, sigma-model and background covariant derivatives are written in terms of the space-time 
and sigma-model metric and of the background fields. 

We use the following conventions. Capital roman indices M — 0, 1, 2, 3, 5 are five-dimensional space-time indexes, 
while greek indexes \i = 0, 1, 2, 3 are restricted to the 4-dimensional Minkowski slices of the space. In this way, we label 
the space-time coordinates as x M = (x 1 - 1 , r), with r the radial (fifth) direction. Lower-case roman indexes a = 1, . . . , n 
refer to the sigma-model (internal) space. 

We write the five-dimensional metric guN with signature — h + + + . We define the five-dimensional connection as 

T P mn = \d PQ ( 9m9qn + d N g QM - d Q g MN ^ , (1) 

and hence the covariant derivatives are of the form 

V M T P N = d M T p N + T P MQ T Q N - T q mn Tq , (2) 

for a (1, l)-tensor, and analogous for other tensors, in such a way as to ensure compatibility with V pgMN = 0. The 
Ricmann tensor, Ricci tensor and Ricci scalars are defined, respectively, as 

RmRN P = ^R^MN ~~ ®M^ P rn + ^MN^ P R.Q ~~ ^ MR^NQ > (3) 

Run = Rmrn R > ( 4 ) 
R = 9 mn Rmn ■ (5) 

One important fact that we will use in this paper is that we will assume the space-time to have one compact 
dimension. It is convenient to choose the coordinates in such a way that the radial direction r± < r < r-i is compact, 
with the slices of space-time with constant r supporting a Minkowski four-dimensional metric. The presence of four- 
dimensional boundaries means that we will need to use the induced four-dimensional analogs of all of the above 
geometric objects, which we will label as g^u, V M and so on. The boundary terms are built starting from the 
ortho-normalized vector N M , defined so that 

g MN N N N M = 1, (6) 
g MN N N = 0, (7) 

which implies that cjmn = Qmn — NmNr. The extrinsic curvature is defined from 

K M n = V M N N , (8) 

as the contraction with the bulk metric 

K = g MN K MN . (9) 

The boundary actions will contain terms proportional to K. 

The field content of the 5-dimension action comprises a set of real scalar fields that we label as <& a with a = 1, . . . , n. 
In a way that is analogous to the space-time metric, we indicate with G a b the sigma-model metric, which ultimately 
encodes the geometric properties of the internal space spanned by the scalars. In the case where the five-dimensional 
system is obtained by consistently truncating some higher-dimensional theory, the structure of the sigma-model is 
determined unambiguously by the fact that the scalars parameterize some coset space, which is in general non- 
compact, and which emerges from the fact that the compactification of the internal space results in the breaking of 
some global symmetry of the underlying theory. We will keep the sigma-model structure as general as possible, hence 
not committing ourselves to any specific realization of this structure. 

Sigma-model indexes are lowered and raised by the sigma-model metric G a b and its inverse G ab defined by G ab Gbc = 
5 a c . When unambiguous, we use lower indices to denote field derivatives with respect to the $ a , so that for example 
given a scalar function V we define 

v a = d a v = dv/d$ a . (io) 
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The sigma model connection is given by 

Q d ab = \c dc ( d a G cb + d b G ca - d c G ab ) , (11) 
and the Riemann tensor with respect to the non-linear sigma-model metric is given by 

na bcd = 9cQ a bd - ddQ\ c + Q a ce G e bd - Q a deQ e bc- ( 12 ) 

Using the sigma-model connection, we define the sigma-model covariant derivative. It is convenient to introduce 
also another convention. When a sigma-model index is placed after a "|", it means that the sigma-model covariant 
derivative with respect to G a b should be taken, which is defined as acting on a (1,1) sigma-model tensor X d a by 

X% = D b X d a = d b X d a + Q d ch X\ Q\ h X A c , (13) 

and analogous expressions for other tensors. 

Finally, one needs to rewrite in a covariant form those objects that have indexes both on the space-time and on 
the sigma model, making them into generalized tensors. However, we do not need to write the general form of 
the covariant derivative for this case, because our theory is written only in terms of space-time scalars (that carry 
sigma-model indexes) and the metric (which is a sigma-model scalar), so that the only object we actually need is the 
background-covariant derivative T>m defined for a five-dimensional scalar, sigma-model (1, 0)-tensor a a , via 

V M a a = d M a a + G a bc d M 3> b a c , (14) 

where i> b means that $ & is evaluated on the classical background (as is Q = £($")). Notice that in the following we 
will assume the background functions to depend only on the radial direction r, and hence only the fifth component 
of the background-covariant derivative has a non-trivial connection contribution, while the other components reduce 
to ordinary derivatives. 



B. ADM formalism 



We derive all the relevant equations using the ADM formalism, the basic idea being that we will rewrite the metric 
9 m n and the scalars $ a as a background function plus general fluctuations, and then decompose the metric in terms 
of four-dimensional tensors by slicing the space-time along the radial direction. We start by writing the metric in the 
form 



9MN + (1 + vf 



(15) 



Because we singled out the radial direction as orthogonal to the boundaries, the normal vector is defined by Nm — 
(0, (1 + v)) and N M = (1 + z/)"^-^, 1), so that g MN = diag{<p", 0} (notice that this tensor is not the inverse of 
9mn)- 

We assume that the background metric satisfies the ansatz 

ds 2 = e^ri^dx^dx" + dr 2 , (16) 

with A — A(r). Similarly for the scalar, we assume that ($ a ) = $ a (r), so that the background depends on the radial 
direction r, but not on x^. We fluctuate the whole system by expanding the scalars as (using the exponential map) 

- CX P4 (^) a EE $ a + ^ - ig^V + ■ • • , (17) 

and the metric (to first order in the fluctuations) as 

9p.v = e 2A {t]^ + hftv), (18) 

with 

h» v = h TT \ + &>e v + d v e" + + ±5» v h, (19) 

where K FTil v is traceless and transverse, and e M is transverse. Altogether, we have the fluctuation variables 
{^,v,v^h TTli v ,h,H,e»}} 



1 Notice that, with some abuse of notation, we identify the fluctuations v and with the components of the metric in the ADM formalism. 
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C. Diffeomorphism invariance 



Here we write explicitly the gauge transformations, and discuss the fact that the boundary actions restrict their 
general form compared to what is allowed in the bulk. The starting point is the five-dimensional diffcomorphisms 



6x M = -£ M , (20) 



which imply 



6<P a = i M d M ^ a , (21) 
SgMN = d M £, R 9RN + d N £ R g M R + £, R d R g M N ■ (22) 

To first-order in the fluctuations, this yields the gauge transformations for all the fluctuations 

<v> = $ ,a e, = drC, w = + e 2A o r e, sh TT t = o, 

Se^ = U^C, SH = 28^, Sh = 6A'C, 

where we defined the projector IP,, = — d , In all of this, all the functions depend on the five coordinates x M . 
The five-dimensional part of the action is going to be invariant under all of these transformations. 

However, we do have four-dimensional boundary actions, the very existence of which is not compatible with all of 
the above. To be more specific, the action is still symmetric under these transformations for generic ^(x^, r), but we 
have to specify how to treat the diffcomorphisms in the fifth direction. Because we will explicitly write the boundaries 
to support localized actions at the points r — ri, with i = 1,2, one must require 2 

rccvo = o. (24) 

This observation plays an important role in the subsequent discussion about gauge-invariance and gauge-fixing. 

D. Gauge-invariant formalism 

Generalizing the notation of [25[ , we define the following variables 

a =V -6A? h ' 

d r (h/A') 

6 ' 

1 (25) 
- 2 9 ^ 



c = e- 2A d^ - 


e- 2A Uh 




6A' 




- d r e^ 







These are a generalization of the Mukanov-Sasaki variable [29(. By inspection, one can verify that these new variables 
are 5d gauge invariant. 

Before proceeding, let us go through a counting exercise. Besides the n scalars in the sigma-model, the off- 
shell degrees of freedom derived from the dimensional reduction of the five-dimensional fluctuations of gMN com- 
prise another 15 components for a total of 15 + n components. Counting in the basis of original fluctuations 
{tp a : v,^,h TTt * u ,h : H : e^} yields the same number n + 1 + 4 + 5 + 1 + 1 + 3= 15 + n, as it should (we counted 
3 for transverse vectors and 4 for generic vectors). The counting of the gauge-invariant variables {a a , b, c, D M , e M „}, 
however, yields n + l + l + 3 + 5= n+10 off-shell components. The five extra components of the original fluctuations 
are pure gauge, corresponding to the diffeomorphisms £ . 

One needs to show explicitly that these gauge-invariant variables are physically equivalent to the original set of 
fluctuations, so that the full set of equations can be rewritten directly in this form, hence removing all the possibly 



2 This expression could be made covariant, if one wanted to manifestly show that one does not need to commit to a specific choice of 
coordinates. 
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spurious gauge artifacts, while retaining all the physical information. In the absence of the boundary action, this is 
straightforward, because one can make a choice of £ M such as to set h = = H = e^, and hence by simple counting 
one can see that the whole system can be rewritten in equivalent form in terms of the gauge invariant variables. This 
can be verified explicitly to hold for the bulk equations (which hence are the ones derived in (25j). 

In our specific case, though, the boundaries restrict the gauge transformations allowed. One has to show that the 
thus constrained system can still be completely rewritten in terms of the variables in Eq. (|25|) . In order to do so, we 
use a different strategy. First, we observe that because of the restriction on £ r in Eq. (j2"4"|) , at the boundaries one can 
define two independent 4d gauge invariant variables, which replace c: 

e~ 2A Dh 

ci - ~ (26) 
c 2 = e- 2A d^ - h) r H , (27) 

so that c = Ci + C2- The original system of fluctuations is certainly equivalent to the set {a, b, Ci, C2, 0, e}. Then, 
we have to show that the boundary conditions actually remove the extra degrees of freedom, hence proving that the 
presence of the boundary actions consistently restricts both the fluctuations and the gauge transformations, so that 
the whole system (bulk and boundaries) can be fully expressed in terms of the variables in Eq. (|25p . 

More precisely, one can show that one of the boundary conditions can be written in a 4d gauge invariant form as 



C2 



0. (28) 



This, together with the fact that in the bulk the gauge freedom allows to always set c = Ci, means that c 2 is actually 
not a physical degree of freedom, and it can be set to zero everywhere, hence allowing for the whole set of fluctuation 
equations and boundary conditions to be written purely in terms of the variables in Eq. (|25[) . In the appendices, the 
boundary condition (1281) is derived by first choosing a gauge transformation £^ such that v^(x^r) = everywhere, 
then showing that at the boundaries d r H\ ri —{). Finally, we can use the residual gauge £ r to set d r H = also in the 
bulk (together with c 2 ). 

Working in the gauge v^{x,r) = 0, and making use of the boundary condition (|2"5j) . there is a straightfor- 
ward and natural one-to-one map between the fluctuations {<p a , v, h, e^, ^i TT ^} and the gauge invariant variables 
{o a , b, c, e M j,}. The main advantage of the gauge-invariant formalism is that it allows to decouple the equations in 
a very simple way pol ]. hence rendering the calculation of the spectrum much easier. In practice, this means that 
the equations for b and c are algebraic equations relating them to the dynamical variables a a , while the equations for 
0^ and decouple. Hence, the spectrum of scalar fluctuations can be identified by solving a set of n second-order 
differential equations involving only the variables a a , subject to boundary conditions that, again, involve only the 
variables a a . 



III. DYNAMICS 



In this section we summarize all the important equations that one has to solve in order to study the spectrum of a 
given model. 



A. The complete action 



We are now ready to write explicitly the complete action. Our starting point is the general definition of the 
five-dimensional sigma-model coupled to gravity. We write the action as 

ii? + £ 5 ($ Q ,dA f $ a ,.9) 
+ ^Kr - ri) [c K K + d^ a ,~g)] 

-^S(r-r 2 )[c K K + C 2 (^ a ,d^ a rg)}} , (29) 

where R is the Ricci scalar, where K is the extrinsic curvature, where the coupling ck = — 1/2 is fixed by consistency 
and where Ci are the sigma-model actions. The step function is defined by O = 0(r — r{) — 0(r — r 2 ). 



d 4 xdr 



-gQ 
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We define the action of the matter fields in terms of the real scalar fields $ a = $ a (ir M , r) as 

1 



£5 



G ab g M "d M $ a d If &-V(& 1 ), 



Ci = -A (1) ($°), 
C 2 = -A (2) ($ Q ). 



(30) 

(31) 
(32) 



Hence, we assume that no kinetic term is present at the boundaries, but only localized potential terms. We will 
provide the explicit forms of the A($) terms later. Notice the different sign with which they enter the complete action. 



B. Background equations 

The background is determined by solving the classical equations, assuming that all the functions defining the 
background depend only on the radial direction r, and not on the x M . We take the variation of the complete action 
in order to determine the equations of motion for A(r) and $(r). 3 The equation of motion for the scalars is 



$" a + 4A'$ /a + g a bc <$>' b & c - V a = 0, 



while Einstein's equations read 



6A' 2 + 3A" = -G a6 $' Q $' b - 2V, 
QA' 2 = Gab*" 1 *"' - 2V. 

The boundary conditions satisfied by the background, at the UV (r = r^) and IR (r = r*i) are 



x a 



G ab d b X (i) 



(33) 



(34) 



(35) 



A' 



- -gA (i) 



(36) 



Thus we see that to quadratic order the localized potentials A(j)(<5°) are constrained to have the following form 4 

1 



X (2) 



3 -A> 
2 



G afc $" 



G afc $" 



($" - + ±D a D b \ (2) (<S> a - $»)($" - , 

r 2 ^ 



(37) 



were <j>" 2 = ^> a | r - 1 2 are the values assumed by the scalars at the boundaries. In effect, this is equivalent to imposing 
the condition that the scalars have fixed boundary values 



1. First-order formalism 
Assuming that there exists a superpotential W such that the potential V can be rewritten as 

V = ±G ab W a W b - ^W 2 , 



(38) 



3 See Appendix [X] for details. An alternative way is to take the variation of the action by linearizing in {h, u, (fi a }. The terms proportional 
to then give us the equations of motion satified by the background in the bulk, whereas the terms proportional to <5 = d r @ give us 
the boundary conditions. 

4 Note that the covariant derivatives of Xu\ arise due to expanding according to the exponential map Eq. I |17l . 
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with W a = d a W — dW/d& a , then the system can be reduced to a set of n + 1 first-order equations 

2. 



A' 



3^ 
= G ab W b = W a . 



(39) 
(40) 



in the sense that all solutions to Eqs. <J39J)- (|40J) are also solutions to the second-order equations Eqs. (|33^ - (j34|) . However, 
for these to provide solutions to the original system, also the boundary conditions must be satisfied, implying a set 
of constraints on the form of the localized potentials A^y The general form of the localized potentials must be 



A ( i) = + W c ($i)(* c - $?) + ^D d D c \ a) {$° - $l)(* d - $f) , 

A (2 ) = W($ 2 ) + W c ($ 2 )($ c -^) + ^D d D c \ (2) ($ c - $ c 2 )($ d - . 



(41) 
(42) 



Before we move onto studying the spectrum of fluctuations, a brief comment is needed. One might, legitimately, 
wonder why we allow ourselves the freedom to add localized potential terms for the scalars, but not localized kinetic 
terms, and/or more general functions not only of the fields but also of their (four-dimensional) derivatives. The fact 
that we truncate A(j) at the quadratic order is just due to the fact that we are interested here only in the equations 
for the background and in the linearized equations for the fluctuations, so that higher-order terms have no effect. On 
the other hand, while terms that depend explicitly on the four-momentum q 2 do not enter the background equations, 
they do enter into the boundary conditions for the fluctuations. However, their presence and structure in entangled 
with the problem of holographic renormalization, in the sense that in order to systematically identify what such terms 
are needed, one has to study the (model-dependent) structure of divergences in the two-point functions. This suggests 
using some caution when discussing the spectrum. 



C. Boundary conditions for the fluctuations 



As discussed in Section III D[ the equations for the fluctuations of the scalars and the metric can be rewritten in 
terms of gauge-invariant fields, up to some subtleties that we will discuss later. Most important, we anticipate here 
that there are two possible ambiguities in this procedure, which is not well-defined for the (interrelated) cases when 
Q' a = and/or when the four-dimensional momentum □ = — K 2 = m 2 = q 2 vanishes. On the one hand, the procedure 
is rigorous even for the case in which these arc infinitesimal, and on the other hand all the examples one might think 
of that have any phenomenological relevance will not contain exactly massless scalars, none of the global symmetries 
being exact. We will come back to these questions later in the next section. 

With all of these caveats, the final result is that the spectrum can be obtained by solving the following second-order 
differential equation for a set of n gauge-invariant scalar fluctuations denoted by a a : 



Vl + 4A'V r 



-2A 



□ 



V] c - n* h J>' h & d 



4(<$> ,a V c + V a $') 16^$' a $' 



3A> 



9A' 2 



0, 



(43) 



with boundary conditions (suppressing the index i of Am) 

-1 (ya 

I 6 ^ 3A' ' ~ ~ 3A 



5\ + e 2A U- 



AV + + < 



4A'$' a - AV$' c1 







) 3 A' 


XV a 6 



„2Aq-1_^_ fya _ AA'fc'a _ q/c\ ( iV ^b 

3A> V l c 



V 3A' 



(44) 



1. Superpotential formalism 
In the special case where there is a superpotential W, we have that Eqs. (|4"3"]l and (|4"4")) become 



8iv r + wa 



W, 



w»w c 
w 



5"e- 2A n 



a c = 0, 



(45) 
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and 



5\ 



A 11 



3 2A a 



1 



2W 



w c w d 
w 



w d w b 
w 



(46) 



respectively. 



D. Boundary masses 

When the superpotential is known, it is convenient to write 



(47) 
(48) 
(49) 

Notice the different sign in the definition. It is convenient to rewrite the bulk equation for the fluctuations in the 
compact form 



N% EE 


w% 


w d w b 
w ' 


a — 
(l)|c = 


w\ 


n + M°c 


a 

(2)|c = 


w\ 


- M°„ 



e' AA (5 a b V r + N\) e AA (S°V r - N b c ) + 5 a c e~ 2A n 



0, 



(50) 



and the boundary conditions as 



5% + e^D- 1 (ml) 
(mi) a b + Nl + 



a W C W b 
c W 



V r a b 

'n 

a W C W d 
c W 



N'l 



(51) 



and 



5 a b -e 2A n-'(mt) 



a W C W b 



W 



V r a L 



{ m l)\ + N%-e* A U-^rnl)\^N\ 



(52) 



respectively. 

The matrices m 2 encode the degree of arbitrariness in the definition of the boundary terms. These matrices do not 
necessarily preserve the same amount of (internal) global symmetries as the bulk sigma-model action. In particular, 
they provide a source of explicit breaking for possible residual internal symmetries of the sigma-model action. Taking 
their entries to be large ensures that possible pseudo-Goldstone bosons (arising by the spontaneous breaking of such 
internal symmetries) are effectively removed from the spectrum. Hence, it is convenient to take the limit in which m| 
are diagonal, and their eigenvalues are all positive and divergent, so that the boundary terms reduce to 



w 



(S d b V r -N d b )a b =*V 



(53) 



One can think of using this procedure also in the case when the physical system does not have explicit physical 
cutoffs. In this case, this procedure can be thought of as a regulator: one computes the spectrum at finite cutoffs, 
and then studies how the physical spectrum changes when extrapolating the values of to their natural boundaries 
(cither towards the UV, or towards a singularity/end-of-space boundary in the IR). 
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IV. APPLICATIONS 



In this section, we apply the formalism developed in the previous sections to a few examples, in order to demonstrate 
how the algorithm that we propose works, as well as to verify that our results agree with the literature when applicable. 
Although the examples we study here are all simple, in the sense that they consist of only one scalar with a canonical 
sigma-modcl metric, the formalism is also applicable to more general cases with several scalars whose sigma-model 
metric is non-trivial. Indeed, most applications, such as various consistent truncations from supergravity, will fall into 
this more general category, and it was with this in mind that the formalism of the previous sections was developed. 



A. Example A: quadratic superpotential and Goldberger-Wise mechanism 

The first example we consider is related to the Goldberger-Wise stabilization mechanism 12], written in the form 
of [l3j]. This example has two advantages. It is peculiarly simple, in that there is only one scalar $ with canonical 
kinetic term, and has a very simple superpotential. And it has been extensively studied in the literature, thus providing 
us with a way to check that the formalism used here yields results that are consistent with other approaches. Also, 
there is a sense in which, in some limit, any physical system near a fixed point resembles it. 

The superpotential is 



(54) 



so that the potential is 



V 



-3 + ^(A 2 



4A)$ 2 



-A 2 $ 4 . 
3 



(55) 



Notice the form of the quadratic term in the potential as a function of the parameter A, which yields the well-known 
result that the five-dimensional mass is related to the scaling dimension of the dual operators by A/|i 2 = A(A — 4), 
and hence provides a natural interpretation for the parameter A. The boundary potentials are 



A(i) 



\(2) 



A 



1 



A*i($ - *i) - - (A- m\) (* - $0 



---_*'-A3 a (3-* a )--(A 



($ _ $ 2 )2 



(56) 
(57) 



The differential equations and boundary terms for the background are 



A' 



A' 
A' 



A n 
1 $ 2 

3 1 



1 



A 



$' + A$i| 
$' + A$ 2 | 



-A$, 
A 



1 ■+ 
0, 

0, 

o, 

0. 



(58) 
(59) 

(60) 

(61) 

(62) 
(63) 



The solution is hence 



'■■2 



$(r) 

- n 

A 



$ie 



-A(i — n) 



1 , $2 

A $1 



(64) 
(65) 

(66) 



where do, and $1 are integration constants. As a result, a big hierarchy between the (physical) UV and IR scales can 
originate from a small value of A (which is protected), and with natural choices of the unprotected $2/^1 ~ 0(1), 
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hence removing possible sources of fine-tuning. The constant clq can always be reabsorbed into a redefinition of the 
four-dimensional units, and hence can be chosen to be ao = 0, so that when $1,2 — > one recovers exactly the standard 
form of the AdS case with unit curvature. It is often convenient to change variable according to r = — In 2, so that 
with Lo < z < L\ one finds that 

ZT - — -(£)*■ «"> 

The value of Lq and L\ are the UV and IR (length) scales, and in this form it is manifest that an exponential hierarchy 
is naturally generated. 



Spectrum 



Before discussing the spectrum, let us make two observations. In the limit in which A — > the fluctuations ip and 
h decouple from each other. Furthermore, in this limit the field $ has trivial bulk dynamics, while the background 
becomes exactly AdS. As a consequence of these observations, the vacuum is determined by two arbitrary, non- 
dynamical quantities $ = $1 = $2 and r% — r\. They correspond to two moduli. In this limit, one expects (at least in 
this semi-classical analysis) the presence of two massless states in the scalar spectrum, associated with these moduli 
parameterizing the space of (inequivalent) vacua. 

The non- vanishing of A yields a background which is not AdS, at least in the IR. Effectively, this corresponds 
to an explicit breaking of scale invariance, and is ultimately responsible for the dynamical stabilization of the finite 
hierarchy ri — r\. As a consequence, one expects the second scalar to stay light even for very large choices of mf, and 
its mass to vanish with A and $1. In the language of the AdS/CFT correspondence, 4 — A is the dimension of a dual 
operator with coupling proportional to $1, the insertion of which breaks explicitly scale invariance. 5 

Because we are interested in the case in which the background is at least approximately AdS, we choose $1 to be 
small. By inspection of the bulk equation, of the background warp factor A, of the superpotential W, and of its first 
and second field derivatives W$ and W$q>, it is apparent that one can expand the fluctuations in powers of 



= A 2 $Je 



2 (T) 2„2Ari 



by writing 



o(r) = ao(r) +xax(r) H , 

and replacing in Eq. (|5P|) . At the leading O(x ), the bulk equation becomes 







(d r - A + 4) (d r + A) + 



-2r 2 

e q 



fl§(r) : 



where the superscript k refers to the heavy KK- modes, and the boundary conditions 

= 



- (mf -d r -A) a k - \m\ < —^- (d r + A) a* 
o q 

9 p 2r W 2 

= - (-ml - d r - A) a h Q + L -mf—^- (d r + A) 4 

o q 



(68) 
(69) 

(70) 

(71) 
(72) 



where we kept explicitly terms proportional to W$ = — A < I>ie~ A ( r '~ ri ) for reasons that will become clear soon. 

The heavy states can be discussed by looking at the solutions to these equations setting W<s,\ Ti = 0. The solution 
to the bulk equations is 



4 = e 2r [dJ 2 -A (e r q) + C2Y2-A (e r q)] , 



(73) 



5 Notice that the unitarity bounds imply that the five-dimensional mass is M^L 2 > —4, which is automatically true provided that the 
superpotential is given by Eq. II 54 I I. This means that when A < 1, the only possible interpretation is that the dual operator has 
dimension 4 — A. 
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with Ci determined by the boundary conditions 



= — (7% — d r — A) do 
= - (-mf, -dr-A) a 



(74) 
(75) 



The details of the spectrum depend on the specific choice of the m 2 terms. Yet, in general the mass gap is related to 
the zeros of the Bessel functions J2-A, at least for r-^ 3> r±. The mass gap is hence ire~ ri = tt/L±, as it is sensible to 
expect. 

Conversely, the lightest states must be treated by keeping the term depending on 1/q 2 in the boundary conditions, 
while dropping the one proportonal to q 2 in the bulk equations. The reason is that in the limit where Wq> — > we 
expect a massless state to be present. Hence, the lightest state will itself have mass 0(x). It is then convenient to 



write q 



xq 2 and W$ 



xe 



-2Ar 



The bulk equation hence becomes, at O(x ) 



(d r - A + 4) (d r + A) a d , 



and the boundary conditions 







„2r(l-A) 

■- > - 2 (d r + A)a d 



i\- d r - A) a; 



(-ml - d r - A) do + 



q 2 

,2r(l-A) 



(d r + A) a d 



(76) 

(77) 
(78) 



where the superscript dp indicates that we interpret this state as a light (pseudo-)dilaton. The bulk equation is 
satisfied by 



c\e 



-Ar 



c 2 e 



-(4-A)r 



(79) 



Notice that, because of the presence of q 2 in a denominator in the boundary conditions, it is sufficient to solve the 
leading-order bulk equation in order to derive the (subleading) O(x) value for the mass of the dilaton. 

One has to solve two algebraic equations to determine the ratio C1/C2 and q. We do so for the extreme case mf — > 00. 
The result is 



_ 2 _ 4e- 2 ( ri+r2 ) (e 2ri - e 2r2 ) (A - 2) 



3 ( e 2(A-2)n _ e 2(A-2)r 2 ) 



which for r\ = yields the dilaton mass 



4A 2 $ 2 : 



A 1 



-2r 2 



3 1 - e 2 ( A - 2 )^ 



(80) 



(81) 



This result, obtained so easily, is in splendid agreement with the literature (23), |24j. 

Notice that in deriving the mass of the light (pseudo-) dilaton we did not make any assumptions about A, aside 
from requiring it to be positive. Hence all of the above holds for generic A, not just for the A <C 1 case. Also, the 
last factor in Eq. (|8Tj) . dependent on r2, assumes unit value for large r-i ^> 0, provided A < 2. It turns negative 
when A > 2, at which point however the negative sign is compensated by the A — 2 factor, ensuring that the mass is 
positive. In this case the mass vanishes for asymptotically large values of Ti. 



2. Zero modes 



We devote this short subsection to analyzing more in detail the limits in which m 2 d 
conditions in Eq. (|77|) and Eq. (|78l) reduce to 



0. For m 2 = 0, the boundary 



(d r + A)a[ 



= 0. 



(82) 



and there is a massless state with bulk profile o|j oc e Ar = z A . In the A — > case this profile becomes constant. 
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The existence of this zero- mode for to 2 — > is a very general property for any system of n scalars. If a superpotential 
description exists, and we apply the formalism of Eq. (f50|). (IBTl) and ([52]) , notice that 



(V 



, W c 

-V r -N a c ) w =0. (83) 



This implies that a a = W a /W always solves Eq. (f50| for q 2 — 0. 

o a satisfies the boundary conditions obtained by setting mf = 0, which reduce precisely to 

( 8 a c V r - iV a c ) a c = 0. (84) 

This observation shows explicitly that there is always a massless state when to 2 = 0, which is the one discussed at 
the beginning of this subsection. Notice however that this is in general the result of fine-tuning of to 2 , and caution 
must be used in interpreting this result. 

In the other extreme, more physical case, in which m 2 — > +oo, the boundary conditions in Eq. (I77p and Eq. (|78[) 
become (for x <C 1) 



3q 

Setting r\ = for simplicity, the solution is 



2e 2r ( 1 - A ) 

(d r + A)a d + a d 



0. (85) 



e 



-Ar /l _ „2(A-l)r 2 > 



< « ^ >r \zr^, > ( 86 ) 

which reduces to Oq oc e _4r + e~ 2r2 in the A — > case (in which this is a massless state). 

Notice a very interesting fact: taking the limit T2 — > +oo in Eq. ()86j) automatically yields a profile that corresponds 
to keeping only the subleading behavior in the generic solution. This shows for a concrete example that the procedure 
we are implementing automatically reproduces the results obtained with the more widely adopted idea of defining 
the spectrum only in the absence of a UV boundary, by imposing that the solutions to the fluctuation equations 
vanish at infinity as fast as possible with r. Finally notice one important fact about Eq. (|86l) . The bulk profile of the 
massless state that is present when A — > is correctly identified by first studying the A ^ limit, in which case the 
corresponding state is light but not massless, and then taking the A — > limit at the end of the calculation. 

3. Discussion 

Let us discuss now what happens for generic values of mf. We keep working under the assumption that x <C 1, 
so that Eqs. (I7T)1) . (|71l) and (jT2"j) . which describe the generic bulk profiles, still hold. Let us consider first the lightest 
state. Solving for generic values of mf, and for simplicity setting r\ = 0, yields 



3 ( e 2(A-2)r 2m 2 ( m 2 + 2A - 4) - {m 2 - 2 A + 4) TO 



2^ ( 87 ) 



2 



This result is completely general (i.e. valid for any to 2 and A). Notice that, as we already know, for to 2 = one 
finds a massless state, and for to,; — > +00 the finite result in Eq. (f5Tj) . suppressed by x. The former choice should 
be avoided, if what one is trying to understand is whether the scenario in question predicts the existence of a light 
scalar, irrespectively of the details of the UV dynamics and of the coupling to other sectors of the complete theory of 
interest (such as the SM). 

Also, a pathology appears when taking \mf \ = 4 — 2 A. The reason for this lies in the way in which we wrote the 
boundary terms. Consider for example the UV term: looking at the coefficient of the ($ — $i) 2 term, replacing such 
pathological choice one has m\ — A = —4 + A. In this case, what is happening is that this is the choice that would 
render massless the excitations around the other solution to the (second-order) bulk equations for the background, 
which has scaling dimension 4 — A rather than A. Therefore, this choice should also be avoided. 

One special comment about unitarity. The way in which we wrote the system, in terms of a quadratic superpotential, 
yields a potential in which the (five-dimensional) mass term Mf = A 2 — 4A > —4 is always above the unitarity bounds, 
irrespective of the value of A. The bound is saturated for M 2 = —4, or A = 2, in which case the theory is close to a 
special transition point that we will discuss at length elsewhere. In proximity of this point, the mass is anomalously 
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suppressed. Yet, we did add boundary terms with arbitrary parameters m|, which distort the spectrum, and one has 
to check that in doing so no tachyon state has been added. One can easily verify that for r-i — > +00, and A < 2 this 
implies that one must enforce the choice to 2 > 0. For T2 — > +00, and A > 2 one conversely must impose m\ > 0. In 
general, given a choice of A and m\ one has to verify that no tachyon is present. 

Finally, we turn to the heavy modes, which have bulk profile a§ up to 0(x) corrections. We already stated (in 
Section IIV A II) that the heavy modes form towers with separation 7re ri . One has to clarify where the second light 
state mentioned several times ends up, for general values of to 2 . The general solution can be obtained by simply 
solving for the integration constants and for q 2 in the equations for o§. This is a somewhat intricate exercise, which 



is not very illuminating. Yet, there exists an interesting limiting case: for mf 



and 



one finds that 



one light state has the mass m? d discussed earlier on in (|8Tj) , while at the same time an exactly massless state also is 
present. Hence, in this case it is clear that two light scalars are present, while the heavy states start appearing with 
masses proportional to 0(ire~ ri ). Effectively, what is happening is that the whole tower is shifted down in this limit, 
and the first excited state in the tower becomes parametrically light. The numerical study performed in the next 
section will make these observations more clear. 



B. Example B: a consistent truncation of type IIB supergravity 

Let us now consider a different example. We still consider the case where only one scalar $ is present, and the 
sigma-model is trivial. But now the superpotential is 



W = — 



so that the potential is 



1 + cosh 2 




V 



-3 cosh 



-3+ -2A 





4 


/a 




H — A sinn 


2\ — $ 


V 3 


8 


V 3 



A? 
2 



(89) 
(90) 



where in the last expression we expanded for small $. Notice how this expansion is in agreement with the potential 
of the previous sections, at leading order. The difference is important only away from the AdS fixed point $ = 0. 
The solution to the bulk equation is 



4> 



rctanh e " Ar+Cl 



(91) 



and we can always choose c\ = for simplicity, setting the radial coordinate in such a way that $ diverges for r — > 0. 
The warp factor is 



A = -L]n(-1 
2A v 



Replacing r = — logz, and expanding for small z — > 0, 



2Ar\ 



(92) 



3> 



(93) 



not surprisingly. Notice however that, as opposed to the what we did in the previous sections, the coefficient in front 
of z A is now fixed. There is no free parameter analogous to $1. The integration constant that we set to zero would 
change this coefficient, but it would also change the position of the singularity in the IR. In this sense, if we think of 
$ 7^ in terms of spontaneous symmetry-breaking, in this model there is a direct link between the formation of such 
a symmetry-breaking condensate and the end of space in the IR (which one would like to associate with confinement). 
For A = 3 the potential becomes 



V = -cosh 2 $(-5 + cosh 2$) 
4 



(94) 



18 



which is the five-dimensional potential obtained by consistently truncating type IIB supergravity on a Sasaki-Einstein 
manifold discussed in [30| in the context of holographic superconductivity, and that has a very long history in the 
context of truncations of type IIB to 5D supergravity intended to yield the duals of controlled deformations of Af = 4 
SYM, being a special case of the GPPZ flows [3jJ|. In this case the full lift to 10-dimensional type-IIB supergravity is 
known. Specifying A = 3, the background scalar is 



and the five-dimensional warp factor is 



<f> = arctanh e ir , (95) 



A = Ilog(-l + e 6r ) , (96) 



where we chose an integration constant in such a way that A — > r far in the UV. Notice how this background is 
practically the same as the one discussed in the previous example, with the choices A = 3, r\ — and $i = 1, aside 
from a very narrow region near the IR boundary, where a singular behavior appears. 

For our purposes, it is interesting to study the spectrum with the background solution given by the first-order 
equations, mainly because of its simplicity, which will help us elucidate on the effect of performing the calculation 
with the explicit boundary terms in presence of singular backgrounds. Unfortunately, because of the singularity in the 
IR, there is no small parameter controlling the VEV to expand in. As a consequence, one must rely on numerically 
solving the fluctuation equations. Yet, we can obtain some semi-quantitative information from our previous results. 
We will perform the numerical calculations at finite ^> r\ > 0, and then extrapolate to the cases where the IR and 
UV cutoffs are removed (ri — > and r 2 — > +oo). 

As long as we choose n 3> 0, the background being hardly different from the previous case, we expect the same 
results to hold. In particular, because the non-trivial departure is localized very close to the IR boundary, we do not 
expect any interesting changes in the spectrum of the heavy modes, for which all the approximations we made should 
still hold, up to a possible overall shift of the spectrum. However, more interesting is the case of the light state. In 
this case, the fact that the VEV of $ diverges near the IR boundary means that the approximations we made might 
not hold. It is hence very interesting to see what happens when r± — > 0. We will study this numerically in the next 
section. 

We conclude with a few comments on the (real) supergravity background generated by this action. It must be 
noted that this system has been extensively studied in the literature, and yields a background that is well-known to 
be a badly singular limit of the GPPZ system, which fails to satisfy even the modest demands of (32j . 

One way of seeing explicitly that a problem is present is the following. The lift to 10 dimensions yields the 
metric [If 6 

ds 2 w = cosh$ (e 2A d^ 3 +dr 2 ) + d0 5 , (97) 

with dQ§ the internal metric, which depends in general on $ and r, besides the coordinates of the five-dimensional 
internal manifold. 

The internal-space structure of the metric is not very important for the present discussion, what matters is that 
we know the warp factor cosh $ needed to lift the five-dimensional metric to ten dimensions. This information allows 
to use the background and compute the Wilson loop. One must solve the classical problem of determining the 
configuration of a (probe) string, the end-points of which are bounded to a D3-brane fixed at some radial position, 
that we can identify with the UV cutoff r 2 [34|. by minimizing the Nambu-Goto action. 

Following the standard procedure (see also [2l|). one first defines the functions f 2 = gtt9xx and g 2 = gttgrr, in 
terms of the elements of the 10-dimensional metric. The separation between the end-points of the string Lqq can be 
computed, by using the auxiliary effective potential V 2 ^(r) = f 2 (r)(f 2 (r) — / 2 (^o))/(s f2 ('")/ 2 ('"o)), as a function of 
the minimum value of the radial direction rg reached by the string in its fall into the radial direction, via 

lqq = 2 r v%- (98) 

We are not going to do this exercise here, but we want to observe the fact that in the case we are discussing one finds 
that 

f 2 = (l-e-^-^-l + e 6 ") 2 / 3 (99) 



The dilaton is constant hence there is no real difference between string frame and Einstein frame. 
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is not monotonic, but rather has a minimum at r = (1/6) ln(3/2). (Also, f 2 diverges at the singularity, which means 
that the model fails to satisfy any of the criteria in |32j |. as anticipated). In turns, this means that the string cannot 
fall all the way into the bulk towards the singularity, but can at most reach down to f. Ultimately, this means that 
the singularity is bad enough that probing the system with extended objects is going to yield unphysical results, and 
hence one should not think of this as a complete model, in which the background captures all the physics of the dual 
confining theory. A resolution of the IR singularity would be needed. 



1. The A = 1 case 

The spectrum of the actual GPPZ system for less pathological cases has been discussed for instance in [33|, and the 
spectrum of this model for A = 1 is discussed for instance in the first reference in (2f|. We briefly digress here and 
redo this last calculation, which can be performed analytically, and which is of marginal relevance to the rest of the 
paper. We apply our procedure, keeping r-i S> t\ fixed, and considering values of r± <C 1, very close to the singularity. 
We limit ourselves to the rrii — > +oo case. By using the fact that N = W$$ — (W$) 2 /W = —A, and specifying A = 1, 
the bulk equation for the fluctuations is 

(-l + e 2r )a'V)+4e 2 V(r) + (g 2 + l + 3e 2r )a(r) = 0, (100) 



subject to the boundary conditions 



4r (d r + 1) a 



(101) 



By solving this equation and imposing the boundary conditions, and then taking the T2 — > +oo limit, and n — > 0, 
one finds that the spectrum is given by 

m 2 = 4n(n + 1) , for n = 1, ■ ■ ■ oo , (102) 

in agreement with the literature. 7 Notice in particular that there are no parametrically light states: the lightest state 
has a mass m 2 = 8, while extrapolating (outside its regime of validity) Eq. (|81[) computed in Example A, with A = 1 
and $i = V3 would yield to 2 , = 4. 

Concluding this short exercise, let us make two comments inspired by the fact that these results agree with the 
literature. First of all, the boundary conditions we use do not rely on the concepts of normalizability and/or regularity. 
They are simply defined algebraically, in terms of the background functions, and there is no need to analyze on a 
model-by-model basis the fluctuations near special points, in order to decide what is physically acceptable and what 
not: the whole procedure is automatically taking care of this, because the boundary actions implement (both on the 
background and on the fluctuations) all the physical requirements. Second, for this particular example, the procedure 
we follow, in which two cutoffs are explicitly present, ultimately yields the same physical results as other procedures, 
once the limit of removing the cutoffs is taken (provided we take m 2 — >• oo). This suggests a general procedure for 
how to calculate the spectrum in the case of when a singularity may be present for the background, assuming the 
IR and UV behaviors of the background are not too pathological, and the limit of removing the cutoffs can be taken 
without difficulties. 



C. Example C: a phenomenological model with cubic superpotential 



Consider now the following superpotential, for one scalar field $ with trivial sigma-model: 

3 A o A 



W = -- - — $ z + — $ J . (103) 
2 2 3$/ K ' 

This amounts to including a cubic correction to the quadratic potential we studied earlier. This superpotential admits 
two fixed points for $ = = and $ = 



7 Note, however, that for generic values of m?, the spectrum would in general be different. 
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The potential, when expanded near the two fixed points, yields respectively 

Vu = -3+^A(A-4)$ 2 +••• , (104) 

2A$? A 2 $f 1 / 4A$?\ o 

Vj = -3-^-^ + -A(A + 4+^j($- < i> / ) 2 + .... (105) 



In order for the $/ to be an attractive IR fixed point, the cosmological constant must be negative and have larger 
absolute value than it has at the UV fixed point. Which is true for A$ 2 > —18, and in particular for any positive 
value of A. 

The solution of the background equations can be obtained by simply integrating 

AW A 

ft.* = — = --*(*,-*). (106) 
We conventionally choose > 0. Then the solution for $ > $/ of $ < is 

$ Q = (107) 

e A(r-r ) _ I ' y * 

where tq is an integration constant. For r > tq the result is a flow away from the UV fixed point $ = 0, towards a 
singularity at r — > ro at which $ — > — oo. For r < ro, describes a flow from asymptotically large values of $ > 
near r — > ro, to the IR fixed point *j when r — > — oo. 

More interesting is the solution obtained when setting the boundary condition so that < $ < */. In this case 

* = e A ( J!!) + 1 > ( ios ) 

where again r* is an integration constant, which in this case has a very different meaning. The flow described from 
* connects the two fixed points, running from $ — > when r > r t , to $ -> when r <C r*. One can easily see that 
$(r*) = $7/2, so that r* is the scale that separates the regimes in which the theory can be described as a deformation 
of the two fixed points, respectively. 

We are interested in studying the spectrum of the theory defined by the background $. Strictly speaking, this is 
continuous. We perform the study by assuming that there exist two hard-wall cutoffs in the IR and UV, such that 
77 <C r* <C rjj . The scales ri t u are clearly spurious, representing cutoffs put in by hand. In particular, the IR scale 
should appear as a consequence of a relevant deformation driving the flow away from the IR fixed-point. This requires 
extending the system, with more scalars being included. Also, the ultimate end-of-space (singularity) will determine 
the spectrum. In particular, it may be that the physics near the singularity, responsible for its resolution, can be 
described only by embedding the model into a full 10-dimensional supergravity, or even superstring theory. This is 
well beyond the scopes of this simple phenomenological model, and hence one should not be too much concerned 
about the dependence of the masses (in particular of the lightest modes) on 77. Yet, it is of interest to perform this 
calculation in order to understand how the spectrum depends on r*, for fixed choices of ri t u- This will offer some 
guidance as to what happens in actual string-theory models describing the field theory RG flow between UV and IR 
fixed points. 



V. NUMERICAL STUDIES 



In this section we present a set of numerical studies of the backgrounds introduced in the previous section. Besides 
allowing us to check explicitly some of the results, this also allows us to understand how good our approximations 
are. 



A. Quadratic superpotential 

We start from the system with quadratic superpotential, considering generic values of A. The first thing we want 
to do is to understand how precise our results in Eqs. ([5Tj) and ([571) are. Figure [T] shows the dependence of the mass of 
the dilaton on the dimension A, computed in the limit m| — > +00. For (not necessarily very) small values of $1, the 
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FIG. 1: Numerical results. Mass m\ of the lightest scalar fluctuation as a function of A. The localized mass terms are divergent 
(m| — > +00). All plots with n = 0, and ri = 2, 3, 5, 10, 15 (with faster fall-off for A > 2 the higher the UV cutoff). Continuous 
lines represent the approximation in the body of the paper, while the points are numerical results. The background is chosen 
to have $1 = 0.2 in the left panel, and $1 = 1 in the right panel. 



agreement of the numerical results with Eq. (|5Tj) is very remarkable (see the left panel in Figure [1} . The agreement 
deteriorates for larger values of $1, yet the qualitative features are preserved. 8 

A set of remarkable physics lessons that can be read directly off Eq. (|8 1 [) are confirmed. 

• There are no tachyons. At least at this level, there is no reason to question the stability of the backgrounds we 
are studying. 

• There is unmistakeable evidence that two very different behaviors appear for A > 2 and A < 2, indicating the 
fact that A = 2 is a very special point of the parameter space, with peculiar physical features. 

• When A < 2, the mass of the light state depends crucially on the dimensionality of the dual operator A, on the 
normalization of the five-dimensional VEV $1 and on the IR scale ri, but not on r 2 , the UV cutoff. With the 
specific choices we made in the plots, it is clear that even at moderate values of r 2 this dependence amounts to 
a subleading effect. 

• The mass to 2 is anomalously light not only for A <C 1 , which is a well-known and studied result [23|, [24| , but 
also when A ~ 2. This is due to the fact that when the limit A — > 2 is taken in Eq. (|8lj) . the result m d ~ 

is suppressed by l/r 2 . 

• For A > 2 the mass is very strongly dependent on r 2 , being exponentially suppressed in the limit of r 2 — > 00 
(taking the limit by holding $1 fixed) with m 2 , ~ |e~ 2 ( A ~ 2 )' r2 (A — 2)A 2 <£> 2 . Thus, provided r 2 is very large, m\ 
vanishes, for all practical purposes, for all A > 2. Yet, this statement is very cutoff dependent, and needs to be 
taken with caution. 

The more general expression in Eq. (|87| taught us a few important subtleties related to these kinds of systems. In 
particular, we already explained that the five-dimensional sigma-model formalism may yield spectra containing many 
light states, which have nothing to do with scale invariance. They might be related to the light techni-pions and/or 
techni-axions of a generic technicolor model, rather than having to do with the techni-dilaton. Of course, because 
all the global symmetries (including scale invariance of the dual theory) are only approximate, the spectrum results 
from non-trivial mixing among all possible scalar bound states. This also implies that one has to be very careful in 
identifying the nature and couplings of the physical (mass eigenstate) states. It is hence a useful exercise to study 
this problem within this very simplified model, by studying explicitly how the boundary terms m 2 distort the scalar 
spectrum. 

This is done in the four panels of Figures [5] and |31 Figure [5] focuses on values of A = 1 and $1 = 0.2, 1,3,5, for 
ri = and r 2 = 5. We keep m 2 — > +00, but vary m 2 . A few very interesting results emerge, which are very general. 



In the right panel of Figure [T] a couple of extra points can be seen for small A. Let us note that these are spurious states, in the sense 
that if one were to take the limit of T2 — > 00, their masses would diverge and they would decouple from the rest of the spectrum. 
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FIG. 2: Numerical results. Mass M of the lightest few scalar states, for A = 1, r\ = 0, r% = 5, plotted as a function of the 
boundary mass rri2, for mi —5- +oo. The four plots differ for the choice of $i = 0.2, 1, 3, 5 (left to right and top to bottom). 



• When TO2 — > 0, one of the masses vanishes. We already explained the reason for this in Section HV Al 

• An interesting level-crossing pattern develops at intermediate values of m2- In particular, this shows explicitly 
how the mixing between the states is very non-trivial. The composition (in terms of the original fluctuations of 
scalars and metric) of the states corresponding to the pseudo-dilaton is in general very complicated. 

• The plots are restricted to the physically acceptable region of parameter space in which m| > 0. However, 
notice how the mass of the lightest state vanishes as a function of mi for m 2 — > 0. If one where to look at 
large negative values of mf, the spectrum would contain a tachyon. This means that in setting up one of these 
models, some attention has to be given not only to how the (super)potential is chosen, to what background 
solutions one studies, but also to which boundary terms are present. 

• It is only for small values of m\ that the spectrum differs significantly from that of the mf — > oo limit. Already 
for m\ ~ 0(1), the spectrum starts to look the same as for mf — > oo, which therefore is a limit that captures 
the generic behavior. Conversely, taking mf to be small should be thought of as a kind of fine-tuning. 

Figure [3] shows another remarkable fact. By varying A, one can see that for small choices of A and of mg, there 
are actually two abnormally light states in the spectrum. By looking at the superpotential, it is immediately evident 
that what is happing in this case is that the two are related to the dilaton and a pseudo-Goldstone boson, the latter 
emerging from the near flatness of the sigma-model scalar potential. As soon as m 2 and A become generic O(l) 
numbers, both these states will in general be heavy, with masses dictated by the general IR scale that controls all 
of the spectrum. However, when A > 2 there is always one light state (whose mass is suppressed by the UV cutoff) 
irrespective of to 2 . 

Let us make a final remark for the reader, who might find it very bizarre that the typical scale of the KK masses 
in the panels in Figures [5] and |3] are so different. We said earlier on that the scale controlling the gaps is simply 
related to n, and hence one might have expected the heavy states to have very similar masses, up to an overall shift 
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FIG. 3: Numerical results. Mass M of the lightest few scalar states, for $1 = 3, n = and r2 = 5, plotted as a function of the 
boundary mass 7772, for mi — > +00. The four plots differ for the choice of A = 0.2, 1.8, 2, 2.5 (left to right and top to bottom). 



controlled by $1. While this is not what the figures seem to show, there is indeed no contradiction, for a subtle 
reason, which ultimately has to do with how we decided to perform the numerical comparison between the three 
different backgrounds, rather than with the physics. The subtlety emerges from the following observation. In all the 
backgrounds we are looking at, at asymptotically high values of r the background is characterized by a vanishing 
small value for $(r), and consequently the warp factor A(r) ~ r + do- We made the choice of setting ao = in 
all cases. This ensures that the backgrounds are asymptotically all the same, with the same normalizations for the 
4-dimensional Minkoski space-time variables. In this way, the value of T2 always corresponds to the same UV cutoff 
scale. However, in the IR the metric is going to change, because $(r) 7^ 0. The larger the value of $1, the larger the 
departure from AdS. And, more importantly, the departure will appear at higher values of r. As a consequence, it is 
not true that by keeping the same value of t\ in backgrounds with different $1 one is introducing an IR cutoff at the 
same scale. In other words, the numerical value of r\ is not an actual physical scale. It can be converted into such a 
scale only as a function of all the other parameters in the background. We will see a much cleaner example of this 
later on in Example C, when talking about backgrounds that describe the flow between two fixed points, for which 
the geometry interpolates between two AdS spaces with different curvature. 



B. Example from consistent truncation 



Here we will study Example B numerically, and focus in particular on the case when the parameter A in the 
superpotential is fixed to be A — 3. As explained before, the model can then be thought of as having a stringy origin. 
The background develops a singularity in the IR at r — 0. However, as shown in Figure 21 the background fields $ 
and A are almost indistinguishable from those used in Example A, with the choices A = 3 and $1 = e~ Ari , with the 
exception of a very narrow region close to the singularity. 

We study how the spectrum depends on where we put the IR cutoff. In particular we focus on what happens in 
the limit of letting r\ — > 0, keeping the UV cutoff fixed. The results are plotted in Figure [SJ where we also show the 
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FIG. 4: The background functions $ (left panel) and A (right panel), used in Example B (black lines), compared to the case 



discussed in Example A, for the choices A = 3 and $i 



-Ari 



(blue lines). 



behavior of the Example A for comparison. There is a tower of KK modes, and a light scalar. Remarkably, despite 
the presence of an IR singularity, all states in the spectrum converge to their respective (finite) values as the IR cutoff 
is taken close to the position of the singularity. Also, while a shift in the masses of the heavy states is clearly visible, 
this effect is very suppressed for the lightest mass. This is a very interesting result, since naively one would expect 
that the lightest state would be the most sensitive to the IR singularity! 
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FIG. 5: Numerical results. Mass M of the lightest scalar fluctuations, keeping the UV cutoff fixed at T2 — 5 while varying the 
IR cutoff n. The right panel shows a detail of the left one, in which the very lightest state is not visible (notice the scale). The 
plots show the comparison between the results of Example B (black +) and of Example A with A = 3 and <5>i = e~ ri (blue). 

In FigureJSJ the dependence of the spectrum on the UV cutoff is shown. As can be seen, the mass of the light scalar 
tends to zero in the limit of infinite T2, whereas the rest of the spectrum stabilizes. 



C. Phenomenological example: cubic superpotential 



This is to be understood as a toy model describing the RG flow between two fixed points. In Figure [7] a set of 
possible backgrounds for this model are shown, with the background functions $ and A explicitly plotted as a function 
of r. The various backgrounds share the same choices of A = 3, $j = 1, r\ = and = 5. The different curves are 
obtained by varying r*, the value of the radial direction at which the nature of <t(r) changes. For higher values of A, 
the kink becomes more localized. 

Notice that for r* —> +oo, the background becomes purely AdS, with curvature radius determined by the IR fixed 
point. Conversely, for r* — > — oo, the AdS background has unit curvature. Notice also that we have chosen the 
integration constant in the warp factor in such a way that (asymptotically) in the UV all the backgrounds become 
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FIG. 6: Numerical results for Example B. Mass M of the lightest scalar fluctuations, keeping the IR cutoff fixed at r\ = 10 
while varying the UV cutoff r%- The right panel shows a detail of the left one, with the line being the appoximation M <x e - ' 




FIG. 7: Sampling of the functions determining the background in Example C, and used in the numerical analysis, plotted 
against the radial coordinate r. Plots obtained with A = 3, <f>j = 1, r\ = and = 5. The backgrounds differ in the choice of 
r* . Notice that the integration constant in the warp factor A is chosen so that the warp factor of all the different backgrounds 
agrees in the far UV (r — > oo). 



identical. This means that the same numerical choice for UV-quantities (such as r%) corresponds to the same physical 
scale. This is not the case for IR-quantities (such as n). We will come back to this comment later on. 

In Figure HI we show the spectrum of masses M for the first few composite states as a function of r* . The spectrum 
has been computed with mf — > +oo. The different figures show the results for different values of A. Focusing on the 
heavy states, we observe the expected behavior: the spectrum of heavy states consists of an infinite tower of evenly 
spaced KK excitations. There is an artificial suppression of the masses as a function of r». This effect was already 
observed earlier on, when discussing the quadratic potential. Since below r* the curvature is not unit, the fact that 
we chose the warp factors of different backgrounds to agree in the far UV means that the same value of r\ yields a 
different physical scale. An alternative way of defining the IR cutoff is to make it be at the point where the warp 
factor A is equal to zero. This then ensures that the IR cutoff is always at the same energy scale, and in this sense 
it is perhaps more natural from a physical point of view. The resulting plots are shown in Figure |H] (for = 1) and 
Figure [TU] (for 3>j = 2). There is nothing particularly deep about this, aside from suggesting that we need to exercise 
some caution when making quantitative statements relating physical scales to each other. 

Much more interesting is what happens at the level of light states. The limits of r* being small or large are associated 
with the physics of the two different fixed points ($ = and $ = <£>/). For small values of r*, the backgrounds are 
very similar to those considered in the case of quadratic superpotential in Sections IIV Al and IV Al and the spectrum 
behaves qualitatvely the same. For instance, in the case of A > 2, there is a light scalar. On the other hand, the 
backgrounds with larger values of r* only deviate away from the fixed point $ = $/ far in the UV, exhibiting walking 
behaviour from the IR cutoff up to the scale set by r*. Therefore, in this case the light scalar has the interpretation 
of being the analog of the dilaton discussed in the context of walking technicolor. This picture is clearest for larger 
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FIG. 8: Numerical results for Example C. Mass M of the lightest scalar fluctuations, for the choices A = 1, 1.5, 2, 2.5, 3, 3 (left 
to right, top to bottom), <E>/ = 1, n = and ri — 5, computed with backgrounds differing in the choice of r*. The last panel 
shows a detail of the left one, with A = 3. 



values of A in which case an interesting crossing structure develops. The crossing structure makes it apparent that 
the nature of the lightest scalar changes radically as r* is varied, being related to either of the two fixed points, for 
small or large values of r* , respectively. 

In Figure [Til the left panel illustrates that for larger values of the light state whose mass is suppressed by the 
length of the walking region (i.e. r*) requires a longer such region to become light. The right panel of Figure [TT1 
shows the UV cutoff dependence of the light states. As expected, the light scalar associated with the fixed point at 
$ = has a mass that is suppressed by the UV cutoff T2 in agreement with the results found for A > 2 in Example 
A. The light scalar associated with the fixed point at $/, however, is unaffected by the value of r2, caring only about 
r* which acts as the UV cutoff for this state. 
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FIG. 9: Numerical results for Example C. Mass M of the lightest scalar fluctuations, for the choices A = 1, 1.5, 2, 2.5, 3, 3 (left 
to right, top to bottom), $>j = 1 and = 5, computed with backgrounds differing in the choice of r», and by choosing n so 
that A(ri) — 0. The last panel shows a detail of the left one, with A = 3. 



VI. FIELD THEORY INTERPRETATION, DISCUSSION AND GENERAL LESSONS 

In this section, we discuss the physical meaning and implications of the examples we discussed in the paper. Before 
we start, we must remind the reader about the two main subjects of the paper. First of all, we are mostly interested in 
understanding under which conditions a strongly-coupled, (quasi) conformal theory admits anomalously light scalars 
in the spectrum. These can be the result of accidental cancellations, but more often are the result of the spontaneous 
breaking of approximate symmetries. Such symmetries can be internal (giving rise to pseudo-Goldstone bosons, such 
as the techni-pions and techni-axions of a technicolor model) or it may happen that scale-invariance is an approximate 
symmetry, in which case the light scalar is a dilaton. 

Second, the framework within which we work is that of gauge-gravity dualities. What we did was to set up a very 
flexible formalism, that allows to study the four-dimensional spectrum obtained from a five-dimensional sigma-model 
of n scalars coupled to gravity, in the presence of UV and IR boundaries in the radial direction. The formalism 
exploits the diffeomorphism invariance of the five-dimensional theory in order to write the (linearized) equations for 
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FIG. 10: Numerical results for Example C. Mass M of the lightest scalar fluctuations, for the choices A = 1, 1.5, 2, 2.5, 3, 3 (left 
to right, top to bottom), $/ = 2 and r2 = 5, computed with backgrounds differing in the choice of r», and by choosing n so 
that A(n) = 0. The last panel shows a detail of the left one, with A = 3. 



the fluctuations and the boundary conditions directly in terms of physical states. This fact allows to reduce the 
complexity of the general problem to a set of n (coupled) equations involving only n scalar fields. 

We applied this formalism to three examples. Up to now, we focused on the technical aspects, showing explicitly, 
on the basis of these three very simple examples, how the calculations are carried out, and what are the main results. 
All the examples involve only one scalar field, with trivial sigma-model metric, and all admit a description in terms 
of a superpotential. The latter is quadratic (in Example A), a simple hyperbolic function (in Example B), or cubic 
(in Example C). 

For particular choices of the boundary potential, all of the examples can be made to contain at least one exactly 
massless mode. Its composition in terms of the fluctuations of the original scalar and gravity degrees of freedom always 
includes a component that couples to the trace of the four-dimensional stress-energy tensor. In other words, there is 
always a field that can be identified with the dilaton. The crucial task is then to identify under which conditions (on 
the bulk dynamics and on the boundary conditions) does this state stay light and keeps (at least at leading-order) 
the appropriate couplings in order to be identified as a physical dilaton. In order to do so effectively, one has also 
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FIG. 11: Numerical results for Example C. The left panel shows the mass M of the lightest scalar fluctuations, for the choice 
A = 3, <&j = 4, and ri — 5, computed with backgrounds differing in the choice of r, , and by choosing r\ so that A(r\) — 0. The 
right panel shows the mass M of the lightest scalar fluctuations, for the choice A = 3 and = 2, computed with backgrounds 
differing in the choice of r*, by choosing ri so that A(r\) — 0, and varying the UV cutoff — 6 (black dots), r-2 = 7 (red x), 
and V2 — 8 (blue +). 



to anticipate the effect that coupling the strongly-interacting sector (dual to the five-dimensional sigma-model) to an 
external weakly-coupled sector has on the spectrum. Hence, one has to make sure that no serious fine-tuning problems 
are present. 

We saw that the boundary conditions are determined by the background solution, up to a certain amount of freedom 
in the choice of the quadratic m? terms. Special choices of the mf may yield very peculiar results for the spectrum. 
In particular, we saw that setting m\ = automatically implies that an exactly massless state is present. However, 
such a special choice is certainly the result of fine-tuning: as soon as coupling to an external sector is added, there 
is no reason to expect that such a choice is preserved. In general, (perturbative) loop corrections coming from this 
external sector are going to yield corrections to m|, which are UV-sensitive (divergent). For this reason, it is most 
interesting to see if a light state exists when mf — > +oo, in which case one can be confident that no fine-tuning is 
hidden in our procedure. In doing so, one is also guaranteed to break any possible approximate global symmetries of 
the sigma-model (and of the dual strongly-interacting theory), so that if a light state exists, it cannot be due to such 
an internal symmetry. However, a word of caution is needed here: if for some physical reason the global symmetry 
of the internal sector happens to be also a global symmetry of the external, weakly-coupled sector of the full theory, 
then one has to treat the matrices mf appropriately. 

Finally, a completely general comment. We said in the introduction that the physics of massive states cannot be 
completely universal, but rather it is necessarily sensitive to mo del- dependent details about the dynamics. If one 
had the exact dual of a well-known and established technicolor model that fits all the data, this observation would 
not matter, one could simply compute the spectrum in the gravity side of the correspondence, and conclude with 
the phenomenological implications of the results. This is unfortunately not the case, in part because no such a thing 
as a standard technicolor model exists, but also for a more subtle reason. Gauge/gravity dualities, in the context 
of dynamical electro-weak symmetry breaking, do provide precise and effective calculation techniques, in the sense 
that the results do not depend on uncalculable 0{\) coefficients, as is the case for four-dimensional estimates with 
strongly-coupled systems. But one faces the limitation that the gravity description is dual to models that are not 
precisely what one wants. In particular, all we are going to say is valid only in the strict large- N limit, and in most 
of the cases some amount of supersymmetry is present. Hence, in spite of the fact that one can get actual numbers, 
rather than order-of-magnitude estimates, for the relations between masses and couplings of the various states, one 
still needs to consider many different models, and try to understand the parametric dependences, rather than the 
actual numbers. 



A. Quadratic superpotential 

We start from Example A, in which case the superpotential is simply quadratic. It must be stressed that all that 
we are going to say is not restricted to this specific model, but rather it will hold also for any other model that is at 
least approximated by a set of controllable perturbations of a conformal theory, provided, in the whole region between 
the two cutoffs that we introduce, the vicinity to a fixed point controls the dynamics completely. We will hence use 



30 



some of the present considerations also in discussing Examples B and C, when appropriate. 

The analysis of the spectrum of the scalar fluctuations has revealed the existence of a number of limits in which 
a parametrically light scalar emerges. Specifically, we found a light scalar if at least one of the conditions A « 1, 
to 2 < 1, or A > 2 holds. The first case has been discussed at length in the literature, and is known to be interpreted 
in terms of a quasi-marginal deformation of the CFT. The second case is the result of a fine-tuned choice, as we saw 
earlier on, and is hence of marginal interest. We hence focus on the last case, A > 2. 

When A > 2, the standard dictionary of the AdS/CFT correspondence implies that an operator O of dimension A 
has developed a VEV. In which case, one clearly expects a massless dilaton to be present. Because of the finite value 
of the UV cutoff T2 that we use, otherwise subleading deformations (such as the insertion of multi-trace operators in 
the dual theory [HJ]) cannot be neglected, and result in a small mass for such a dilaton. The irrelevant nature of such 
deformations explains the suppression of the mass m 2 , in Eq. (|8ip as a function of r2, and the fact that a massless 
state can be recovered in the r2 — > +00 limit. 

In order to make this more quantitative, let us assume that the dual theory is such that dimO 2 = 2dimO (e. g. 
in the large- N limit). In the study performed in the body of the paper, we take the limit r-i — > +00 by holding the 
VEV $1 fixed. A simple toy-model description of the dual effective potential, due to the multi-trace deformations 
that break conformal invariance for finite r-i, can be approximately given by A 4 ~ 2A (O-vf, where v is the VEV in 
the field theory, and the scale A is related to the UV cutoff. For A > 2, this effective potential is suppressed by the 
UV scale, because the operator O 2 is irrelevant (together with any even higher order O n correction), hence explaining 
our result that m 2 , oc e -2(A-2)r 2 from Eq (Jgjj 

B. Example B 

The second example we considered can be thought of as a completion, in the context of Type-IIB, of the case with 
quadratic potential in Example A, for A = 3 (and $1 = 1) and for A = 1 (and <!>i = y/3). This is nice for several 
reasons. First of all, it means that the superpotential is known beyond the quadratic level, and hence can be used 
with some degree of confidence also far away from the UV fixed point. Second, if one chooses dTl^ in such a way that 
for <J> = it yields the metric on the five-dimensional sphere, one has a precise mapping in terms of deformations of 
Af = 4 super- YM. In the A = 3 case, we are enforcing a non-trivial gaugino condensate [3lj], while for A = 1 we are 
giving a mass to the fermions. Yet, the presence of a singularity in the IR of the background means that this is not 
the complete dual to the four-dimensional field theory, but that some ingredient is missing, so that the physics of the 
four-dimensional theory is well-captured only away from the singularity. In practice, this means that the truncated 
system yielding our very simple superpotential is actually too simple: it does not capture some non-trivial properties 
of the strong dynamics of Af = 4, taking place very far from the UV fixed-point. This is all well known, and discussed 
in the literature. 

It is therefore interesting to understand what happens when computing the spectrum using our algorithmic proce- 
dure, which implies adding cutoffs both in the UV and in the IR. There is not much to say about the A = 1 case. 
This just provides a very nice cross-check, showing how the regulator procedure we use yields the correct results, in 
a simple case in which other arguments can be used to discuss the spectrum without introducing any regulators. 

We hence focus here on the A = 3 case, which we extensively studied numerically. Since A > 2, for finite UV cutoff 
the interpretation is that we add an irrelevant deformation to a theory and enforce a VEV. Provided the IR cutoff r\ 
is far away from the singularity, the VEV is spontaneously breaking scale invariance, as in Example A. However, as 
n is chosen to be close to the singularity, there is no obvious sense in which the theory is still close to a conformal 
fixed-point, and hence no a priori reason to expect the light state to persist. It is hence very interesting, and somewhat 
surprising, that it does. Notice by contrast that the spectrum of heavy states is shifted, to testify of the fact that the 
singularity is not a negligible correction to Example A. Even more, remember that the study of the Wilson loops we 
briefly sketched shows a very dramatic effect, to the point that the string probe cannot even approach the region in 
the immediate proximity of the singularity. 

The fact that the IR is badly singular clearly signals that the strong dynamics of Af = 4 involves other non-trivial 
effects, which are not captured by this simplest truncated model, as we said earlier. Yet, the boundary action localized 
at r = n acts as an IR regulator, which effectively removes from the calculation of the scalar spectrum the pathological 
effects of the IR singularity. The resulting spectrum has all the sensible features expected in a healthy field theory. 
It would be very nice to know if (and to what degree) the spectrum we computed is in quantitative agreement with 
what is obtained in a modification of the model such that the IR singularity is resolved. In particular, we found the 
presence of a very light state (which did not exist for A — I), and one should test whether this state is still present 
in the more complete analysis, rather than being an artifact due to the combined effects of the bad singularity and 
the IR regulator. 
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C. Cubic superpotential 

The field theory motivation for studying this example is that one might be interested in studying four-dimensional 
models in which the RG-flow of a confining (UV-complete) theory happens to come very close to a non-trivial fixed 
point at intermediate energies. One such example of phenomenological relevance is the class of models that embed 
walking technicolor into extended technicolor, in order to explain in a unified picture the three (superficially conflicting) 
requirements of generating large masses for standard-model fermions while at the same time suppressing new-physics 
contributions to FCNC processes and preserving the UV-completeness of the theory. 

There is no known example of a dynamical model the exact gravity dual of which has all the features required 
in the walking technicolor framework. There exist models that reproduce the flow of a confining theory, and there 
are models that describe the flow between two fixed points, but there are no models in which the IR fixed point is 
only approximate, and ultimately the theory confines. Also, models yielding the flow between two fixed points have 
a tendency to be very complicated, such as the Pilch- Warner dual of the flow from J\f = 4 to the Leigh-Strassler fixed 
point, for which the actual background is known only numerically (and within this framework, flows that approach 
the IR fixed-point but do not reach it are described by backgrounds which are badly singular). 

Yet, it is of general interest to know whether, in such a theory, it is at all possible that a light dilaton is present. 
Notice that the answer to this question is not obvious: even when considering the exact flows between two exact fixed 
points, (in which case both the far-UV and deep-IR effective descriptions are provided by CFTs) the theory as a whole 
is not scale-invariant: there exists a physical scale (connected to the p* in our study) that separates the regimes in 
which either of the two CFTs provides a sensible approximation for the physics. So, the actual mass of the lightest 
scalar (the would-be dilaton) will in general depend on this scale, on the two CFTs living at the fixed points (the 
spectrum of dimensions) , and on the specific properties of the flow. 

The study we performed shows that indeed there is a light scalar, under quite general conditions, and that its mass 
depends on the dimension A and on p* . Let us start from the dependence on the dimensions. At the IR fixed point, 
all the (active) scalars must correspond to irrelevant deformations (otherwise the flow could not reach such a fixed 
point). The only important distinction comes from the value of the dimension A, which characterizes the flow near 
the UV fixed point. 

What is most remarkable, is that when A > 2 the dependence on r* of the mass of the lightest state is not a 
monotonic function. There exists an actual maximum of the mass, obtained for non-extreme values of r* . The reason 
for this is that when r* is so large, or so small, that the theory is effectively always very close to one of the two 
fixed points, the mass of the lightest state is suppressed parametrically, because the theory is effectively very close 
to conformal. For intermediate values of r* , the theory is not well described by a small departure from a CFT, and 
as a result a non negligible mass is generated for the lightest state. However, the r* dependence of such mass must 
interpolate between the two extremal values of r* (the UV and IR cutoffs), near which the mass practically vanishes. 
Hence, there is a maximum for such mass. It is curious to notice that numerically such a maximum is (at least in our 
examples) still significantly suppressed in respect to the typical scale of the heavy states, although the actual physical 
relevance of such a fact is questionable. 

VII. CONCLUSIONS AND OUTLOOK 

We conclude the paper by summarizing our main results, by critically reviewing the limitations of this approach 
and by suggesting a set of physically interesting applications. 

A. The algorithm 

Let us first summarize the algorithm to be used to compute the spectrum. We start with a five-dimensional 
sigma-model consisting of n scalars coupled to gravity. For simplicity, suppose that a superpotential W is known, 
and furthermore that there is no obstruction to taking mf oo. The spectrum of scalar bound states can then be 
computed by applying the following steps. 

• Write the action in the form of Eq. (f2"5)) . with the bulk dynamics of the n scalars in the form of Eq. (J3UJ) , and 
the background metric in the form of Eq. (fH)|) . 

• Introduce a UV and an IR cutoff, by writing boundary actions as in Eq. (|3~Tj) and Eq. (132 j) . 

• Determine the background, by solving Eq. f|39[) and Eq. (|40l) . subject to the boundary conditions in Eq. |35|) 
and Eq. ||3BJ>. 
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• Obtain the spectrum by solving Eq. (|5P|) . and then identify which values of q 2 = □ allow to satisfy the boundary 
conditions in Eq. (|53|) . 

• If possible, and physically meaningful, take the limits r2 — > +00 and r% — > ro, where ro is the end-of-space in 
the IR. 

If the superpotential is not known, or it is not legitimate to take m 2 — » 00, all the necessary changes to be 
implemented in this procedure are explained in the body of the paper. It is also implicit that one should familiarize 
oneself with the notation, which is explained in detail in Section UH where all the relevant information is provided 
explicitly. 

B. Limitations of the algorithm 

The algorithm we identified fails, or needs to be partially extended, in the following cases. 

• The bulk action contains terms with more than two derivatives of the scalars and/or the metric. In this case, 
the whole procedure has to be rethought from scratch. 

• The dependence of the two-point functions on the UV cutoff requires introducing q 2 -dependent boundary terms 
in the UV, in order to regularize the theory and remove the UV cutoff itself. If these terms are polynomial, they 
just result in a comparatively harmless modification of Eq. (|53p . to be dealt with via holographic renormalization. 
If they are non-polynomial, then the whole physical meaning of the spectrum becomes questionable, and the 
best thing one can do is to consider the dual theory as some phenomenological model with a physical cutoff T2 , 
which cannot be removed. 

• There are exactly flat directions in the supergravity potential, connected with moduli of the field theory. In 
this case, there are exactly massless states which have nothing to do with the dilaton, and at the technical level 
the D _1 operator appearing in many equations is badly defined. One should find a (mo del- dependent) way to 
overcome this difficulty, either by adding a perturbation that lifts the flatness of the potential, or by further 
truncating the sigma-model in such a way as to decouple the potentially problematic (inactive) fluctuations. 

• In the IR, the space ends in a naked singularity at ro, and the geometry near the singularity is so bad that 
(super)gravity cannot be trusted. In this case, one has to keep the IR cutoff r\ > r , and firm physical conclusions 
cannot be drawn in full generality, until a resolution of the singularity within a more general sigma-model is 
found. However, note that in the one example that we studied where a naked singularity is present, i.e. Example 
B, our algorithm in fact yields finite results. 

• The IR is not singular, however the end of space is known to be described by extended objects that go beyond 
the (super)gravity approximations. Again, in this case one can only keep r± > ro, using a cutoff that chops off 
the space at a scale where the (super)gravity description still holds. It is not known (and it would be interesting 
to know) whether a relation between the spectrum computed with the present algorithmic procedure and the 
actual physical spectrum computed by fluctuating the whole background exists. 

C. Physics lessons 

The main physical motivation of this work, as explained in the introduction, is to understand what kind of (confining) 
strongly-coupled theories yield potentially light scalars in the spectrum, one (linear-combination) of which has the 
couplings of a four-dimensional dilaton. None of the examples we considered is the exact dual of such a theory, and 
yet in all the examples a set of very general lessons appears to emerge coherently. We summarize here these results. 

When the dynamics of the theory is well described in terms of the properties of the RG fixed points and their 
proximity, a light scalar appears to be present, irrespectively of model-dependent details, provided the dimension A 
(defined in the body of the paper) satisfies either A > 2 or A < 1. In the latter case, this is the well-known fact that 
a very small A means that a (quasi-)marginal operator is deforming the CFT, hence introducing a parametrically 
small explicit breaking of conformal symmetry. This is true even in the Higgs sector of the Standard Model, provided 
all the couplings are very weak, and is hence not particularly interesting. A lesser known fact is that this is true 
also when A > 2. Provided factorization holds, one can interpret this particular result in terms of the spontaneous 
breaking of scale invariance via the VEV of an operator O of dimension A > 2. The fact that the (double-trace) 
deformation O 2 is irrelevant ensures that the light scalar present has a mass suppressed by the UV cutoff. 
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However, one has to exert some caution in using these results. First of all, they hold only in the limit where 
factorization is exact, and no non-trivial operator mixing is present. Most important, they depend on what is meant 
by the dynamics being well described in terms of its properties in proximity of the fixed-points. For example, what 
about QCD (or, maybe better, Yang-Mills)? It is well known that QCD admits a trivial UV fixed point, and that the 
departure from it is induced by a quasi-marginal deformation. And yet, there is no light scalar state in the glueball 
spectrum. The reason is that the RG flow goes very far away from the fixed point, and at the confinement scale there 
is no sense in which the theory is still close to it. Indirectly, this also means that it is very unlikely that the physics 
near the confinement scale is controlled by a non-pcrturbativc fixed point: in this case, as we saw, in spite of the fact 
that the whole theory is not scale- invariant, a light state would still be expected, which is not the case for QCD/YM. 

A set of minor caveats emerged during the technical calculations. Summarizing their implications: in the cases 
where we identified the presence of a light state in the spectrum, one should not immediately conclude that such a 
light state is there in general, but should then ask whether the process of coupling the strongly coupled sector to 
external observable sectors still preserves this result. The whole problem of holographic renormalization is implied. 

Finally, we conclude with some important comments about walking technicolor and similar theories. In this case, 
one is interested in a model where the RG flow starts near a fixed point, evolves towards a different fixed-point, but 
then never reaches the latter, and ultimately confines, with the formation of non-trivial condensates that break the 
global symmetries of the theory. What did we learn about this scenario? While we do not have the gravity dual of a 
walking theory, and hence a firm conclusion should be postponed to a future study, yet by combining the results of 
Examples A, B and C we obtained a set of important indications of the presence of a very light state in some cases 
of walking dynamics. 

Let us summarize here all the elements leading to this conclusion. In studying Example C, we considered the case 
where the theory is well-approximated by the dynamics near the UV fixed point for r » r*, and by the dynamics 
near the IR fixed-point for r -C r*. We chopped off the space in the IR, assuming this to be the end-of-space of 
the model. When r* is small, effectively the theory is never really well described by the IR fixed point, and for all 
practical purposes, this is not a walking theory. When is large enough that a sizable range < r < of the 
fifth-dimension can be thought of as describing walking, something interesting happens. There is always a light state, 
that has a natural interpretation in terms of a dilaton. The mass of such state is suppressed by r»: the larger r* is, 
the lighter the dilaton. Ultimately, this result agrees with the expectations from Example A, in the sense that if one 
expands the background near the IR fixed point, the resulting A is always going to be large. From the field theory 
point of view, this must be the case: if the RG flow was governed by a relevant coupling, it could not approach the 
fixed point, which implies that by expanding near the IR fixed point one must find a large value for A. As a result, 
r* represents the scale of explicit breaking of scale invariance, and enters the effective description near the IR fixed 
point by suppressing the coefficient of the irrelevant operator O 2 inducing the explicit breaking itself. This is very 
interesting for practical purposes: it means that one can write the mass of the light scalar in terms of r* and A > 2, 
hence providing a functional relation between the walking scale associated with r* and the mass m 2 , oc e - 2 (A-2)r» 

All of this means that the very fact that the theory walks (in the sense of coming very close to a IR fixed point) 
implies that a potentially light dilaton is present, and that the UV scale at which walking ends effectively suppresses 
its mass. One has then to ask what is the effect on the mass due to the fact that in the IR the theory confines. 
Example B yields some very interesting information, by looking at the A = 3 and A = 1 cases, and comparing the 
results with Example A. Even when the IR ends in a singularity, the mass of the lightest scalar is in substantial 
agreement with what was found in Example A (notice that this is not the case for the heavy KK state, for which 
the shift in mass is significant). The most important parameter appears to be A. If A > 2, effectively the confining 
behavior can be attributed to the formation of a condensate that breaks spontaneously scale invariance, with explicit 
breaking accounting only for a small effect. On the contrary, for A < 2 no light state exists (unless A <C 1 which 
we already commented about) , and scale invariance is broken explicitly. In a realistic walking technicolor theory, the 
condensate that takes the theory away from the IR fixed point is presumably related to the electro-weak scale. On 
the basis of phenomenological considerations it is usually believed that the dimension of the chiral-symmetry breaking 
condensate must be 2 < A < 3, which would imply the existence of a light dilaton in the spectrum. 

This opens a possibility that in certain walking technicolor theories, in which the walking scale is parametrically 
higher than the electroweak scale, a parametrically light scalar is present in the spectrum. This light scalar should 
then be interpreted as a light dilaton, and hence have couplings that are very similar to those of the SM light Higgs. 
Hence, the discovery of such a light scalar at the LHC might be interpreted as a first indication in favor a walking 
technicolor origin for electroweak symmetry breaking. Only the (non)detection of somewhat heavier resonances, in 
the TeV region, would allow to solve this possible ambiguity in experimental signals shared by the weakly-coupled 
Higgs models and the strongly coupled walking models. It would be very nice to test this conclusion for an explicit, 
complete model of the gravity dual of a walking theory, hence putting this line of arguments of firm grounds. 
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D. Outlook 



The whole machinery we put in place and summarized here has its natural application in the study of the spectrum of 
a confining theory the full dynamics of which is well-captured within a five-dimensional sigma-model which is obtained 
by consistently truncating a fundamental, UV-complete theory. No phenomenologically useful such constructions 
exists yet, and hence an obvious direction for further research is to try to identify such a model. Doing so, and 
hence exploiting in full the potential of gauge-gravity dualities, has the great technical advantage that very non-trivial 
properties of the dual theory can be computed, and the algorithmic procedure we provided renders the extraction of 
the spectrum a relatively harmless technical matter. 

In particular, our comments on walking dynamics and walking technicolor are not supported by very robust argu- 
ments, but rather based on circumstantial evidence emerging from a set of toy models, which we chose mostly on the 
basis of their simplicity. Yet, our work shows that it is important to find complete duals of walking theories within 
superstring theory, and that the study of the spectrum of such models could yield very non-trivial results of utmost 
phenomenological and theoretical importance. 

From a more formal point of view, and in the short-term, a large number of interesting questions have been left 
open by this study, and require applying the procedure we outlined to less ambitious five-dimensional backgrounds. 
We mostly concentrated here on semi-realistic toy-models. It would be very interesting to apply all of the above to 
the many well-known backgrounds that exist in the literature, and whose dual four-dimensional theories are well- 
understood. One could then compare the results obtained using the IR and UV regulators to the results obtained 
with other techniques, both on the field-theory side and on the gravity side of the correspondence. This would yield 
important tests of the correctness of this procedure, and might help to shed light on the physical implications for 
models where the field theory is understood only in part. Also, it would be interesting to implement all of the above 
within the systematic program of holographic renormalization. 
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Appendix A: Equations of Motion in the ADM formalism 



Here, and in the following appendices, we derive all the equations of motion, and the boundary conditions for the 
background as well as the fluctuations. Much of the notation is the same as that of [251 ] . 

Starting with the action of the non-linear sigma model with boundary terms (where in the body of the paper d = 4) 



S = / d d xdr 
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we derive the equation of motion for the scalars 
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and Einstein's equations 
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where g^ = g^, and g rt i = g r r = (where the indices /i and v run over the d-dimensional space-time), and tilde is 
used to refer to d-dimensional quantities. 
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We will now rewrite these equations of motion using the ADM formalism. We start by writing the metric on the 
form 



9mn 



n u n^n^ + n 



» + n 2 h 



(A4) 



where comparing to the notation used in the body of the paper, we have that n = 1 + v and — v^. The inverse 
metric is given by 



1 / n 2 g^ + n^n v -n v 
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The tangent vectors are given by X 1 ^ = and X^ = 8". We have a normal vector Nm = (0, n), N 
n _1 (— n^, 1). The second fundamental form is 
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(A6) 



where we use the notation of [25J (which only differs slightly from the one used in the body of this paper i.e. 
Kmn = VjvfiVjv), so that expressions can be easily compared. One can derive the following relations 
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Finally, we are ready to write down the expressions for the equations of motion using the quantities defined above. 
The equation of motion for the scalars becomes 

\d 2 r - In^d^dr + n 2 V 2 + n"n v V„d u - (n££ + d r Inn - n^d^ lxin)d r + 
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Einstein's equations separate into normal, mixed, and tangential components, obtained by projecting with p MN 
n m n n _ g^ X M X N^ pMN = N m x n^ and pMN = xjfX?, respectively. The normal component reads 
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In deriving this expression, the following relations are useful. We have that = ^^-y g M!y , P r/J = — t - , and 

P rr = so that P MN g MN = 1 - d and P MN g M N = - d. Furthermore, g MN g MN = ^ + d. The mixed 

component is given by 
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and the tangential component is 
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At zeroth order in the fluctuations, the equations of motion for the scalars are 

$" a + dA'$' a + £ a bc $' b $' c - V a = X (i) S h 



and Einstein's equations yield 



d(l - d)A' z + 2G ab ^' a <P' b - AV =0 



A" + dA' 2 



d-1 



V = — 



Writing the background as 



we obtain the boundary conditions 



$'(r) = (6(r - n) - 6(r - r 2 ))$'(r)> 
A'(r) = (0(r - n) - 9(r - r 2 ))i'(r), 



.4' 



A (i) 



where we have dropped the hats. 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 



Appendix B: Linearized Equations of Motion 

Let us now expand the equations of motion in fluctuations of the metric and the scalar fields to linear order. As 



explained in Section III Dt we can work in a gauge where 



0. Furthermore, since we want to consider only 



spin-0 fluctuations, we can put e M = h TT> ^ — 0. This leaves us with the fluctuation variables {ip,v,h,H). Thus, we 
expand equations (|A8j) . (|A9[) . (jAlOj) . and (jAllj) using the following rules: 



n = 1 + v, 



In doing so, we will make use of the relations 



n££ = -A' 6j} - 



m = — 



d-i 



2(d- 1) 
-~ 2A Uh 



□ 



2 □ 
d-2 



2{d-l) " 2{d-l) 

R= -a~ 2A Uh, 

which are true to first order in the fluctuations. The linearized equation of motion for the scalar becomes 



d 2 ip 

<i' a 



*U<p a + dA'd r <p a + 2G a bc $' b d r <p c + d d g a bc $ ,b i> V 



dV a 



d<p< 



^-f)d r h + \d r H - d r ^j - 2v a v = s i( x h u + d * x h^- 



(Bl) 



(B2) 



(B3) 
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At first order in the fluctuations, the normal component of Einstein's equations gives 

4& a (V r cp a ) - W a if a - dA'd r h - (d - l)d r H - Wv - e" 2A Uh = 0, 
whereas the mixed component gives 



(B4) 



(d - l)A'v - -d r h - 2$' a Lp a = 0. 
From the tangential component of Einstein's equations, we obtain 



-2A 



rt 2 h rlA' A' 

° r r + —d r h + —d r H - A'd r P + 

2{d-l) d- 1 2 r 2(d-l 



■Dh- 



d-l d-1 



— ^2S l (\ {l) v + d a X {i) tp a ), 



(B5) 



(B6) 



and 



Equations (|B3|) . (|B6|) . and (IB6|) lead to the boundary conditions 



(B7) 



(B8) 



d r h 



2(d — 1) 



and 



a.iJ = 0. 



(B9) 



(BIO) 



Equation (|B8[) gives the boundary condition for the scalar fluctuations. In the special case of one scalar, this expression 
agrees with the one given in (24|. (|B9[) is actually implied by equation (|B5[) and therefore does not give us any new 
information. Finally, the boundary condition for the variable H shows that we may put 



C2 



-1A 



d ll u"--d r H 



= 



(Bll) 



at the boundary. The reason is that since the form of the boundary conditions must obey 4d gauge invariance, 
the expression H' = must generalize to C2 = had we included the ^ M fluctuations as well (this can be checked 
explicitly) . 



Appendix C: Translation to Gauge Invariant Variables 



The significance of that one of the boundary conditions (IB11[) reads C2 = is that using this relation, and the 
fact that C2 can be gauged away in the bulk, one now finds a one-to-one map between the gauge-invariant variables 
(a a , b, c) (defined in (1251) ) and the fluctuations (ip a , v, h). We have that 



h = -2(d- l)A'e 2A n _1 c, 
tp a = a a - $'°e 2j4 n- 1 c, 
v = b-e 2A D- 1 (2A , c + d r c) 



(CI) 



Let us proceed to derive expressions for the equations of motion in the bulk and the boundary conditions in terms of 
(a a , b, c). The boundary condition for the scalar fluctuations (|B8I) becomes 



V r a a 



= $' a b + ( V a - dA%' a - AV$' C ) e M n _1 c + A a , c a c 



(C2) 
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For the equation of motion for the scalar fluctuations in the bulk (|B3[) . we obtain 



Vi + dA'V r + e~ 2A D 



ya _ a fc'bfald 

v c bed 



a c -& a (c + d r b)-2V a b = 0, 



whereas the linearized Einstein's equations (|B4[) . (|B5I) . and (IB7|) lead to 

2 / = ,_ 4F#'„ 



(C3) 



(C4) 



and 



(d-l)A" 



S r c + dA'c - e- 2A Db = 0, 



(C5) 



(C6) 



resp ectively (where we have used the latter two equations in deriving the first). We recognize these equations from 
[25( | and [23]. Solving for b and c in (|C4|) and (|C5[) . we can now rewrite (|C2|) and (|C3j) in terms of only the scalar 
fluctuations a a . We obtain that a a satisfies the following equation of motion in the bulk 



V 2 r + dA'V r + er 2A U 
with boundary conditions 



a , 4(<f ,a V c + V a V c ) , \WV a V c 



v\ c - n a bcd & b <s> 



(d-l)A' {d-l) 2 A' 2 



a c = 0, 



2$1 



X" 



2$' a $^ 



„2A 



LT 



"I ? 

d-lU' V 



4V-K 



!b (d-i)A' 1 " ~ (d-i)A'V i c ^yv(d-i)^' 

In the special case where there is a superpotential W, these expressions can be written as 

W b W r 



Vb 



s%v r + w? h — ^ - ^^wst 1 ) ( s b c v r - wh 



w 



d-i 



w 



8 a n e- 2A U 



and 



5\ + e 2A U~ l ( A a , - W a ,. 



W c W b 



w 



v 



w 



1 (n 



-wr 



w c w d 
w 



w d w b 
w 



(C7) 



(C8) 



(C9) 
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